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RESUMEN 

En los edificios abiertos el centro de rigidez esta bastante distante del centro de masas 
por lo que durante un sismo tienen problemas de torsion. En efecto los porticos debiles se 
desplazan mucho mas que el centro de masa y en algunos casos los porticos fuertes tambien 
se desplazan mas que el centro de masa. Todo esto ocasiona que los edificios con estas 
caracterfsticas tengan un mal comportamiento sfsmico, si se construyen en la forma clasica. 
Ahora si estos edificios se construyen sobre aisladores de base se reducen notablemente los 
problemas de torsion ya que los desplazamientos laterales en la superestructura son bastante 
bajos. 



En este artfculo se analizan 20 edificios de hormigon armado de 3 a 6 pisos, con 
diferentes dimensiones en planta. Se los analiza ante la accion de 28 sismos artificiales que 
son compatibles con los espectros del Codigo Ecuatoriano de la Construccion para los cuatro 
tipos de suelo estipulados en el Codigo. Estos edificios se analizan con y sin aisladores de 
base, se encuentran los factores de amplificacion por torsion en los porticos extremos y se 
propone una nueva forma de cuantificar estos factores en la que se refleja la bondad del uso de 
aisladores de base elastomericos en las construcciones. 



ABSTRACT 



In the open buildings the rigidity center is quite distant from the center of mass being 
the reason why during an earthquake they have torsion problems. In effect the weak porches 
move much more than the center of mass and in some cases the strong porches also move 
more than the mass center. All this causes that the buildings with these characteristics have a 
bad seismic behavior, if they are constructed in the classic form. Now, if these buildings are 
constructed upon base isolators, they reduce the torsion problems remarkably since the lateral 
displacements in the superstructure are quite low. 



This article examines 20, 3 to 6 buildings of reinforced concrete, with different 
dimensions on site. They are analyzed under the action of 28 artificial earthquakes that are 
compatible with the spectra of the Ecuadorian Code of Construction for the four types of land 
stipulated in the Code. These buildings are analyzed with and without base isolators, the factors 
of amplification by torsion in the extreme porches are found, and a new form is proposed to 
quantify these factors in which the kindness of the use of elastomeric insulator is reflected in the 
constructions. 
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1. INTRODUCCION 



Un edificio abierto, es un edificio esquinero cuyas paredes que colindan con las 
estructuras adyacentes, estan acopladas a la estructura. En la figura 1 se muestra un edificio 
abierto de dos pisos, los porticos que no tienen mampostena dan a la calle y los porticos con 
mampostena han sido construidos incorporandoles al edificio con chicotes o se construyo 
primero las paredes y luego las columnas con lo que la mampostena queda acoplada a la 
estructura. 




Figura 1 Descripcion de un edificio abierto 



Normalmente los edificios abiertos son utilizados para locales comerciales y la primera 
planta tiene una altura mayor a las otras plantas. En la figura 2 se muestra un edificio abierto 
ubicado en la ciudad de Armenia, que colapso durante el sismo que afecto al eje cafetero de 
Colombia el 25 de enero de 1999 y que tuvo una magnitud M b =5.9. Aguiar (1999). En la 

figura 2 se aprecia que el tercer piso colapso completamente sus columnas no resistieron las 
cargas debidas al sismo. 





Edificios abiertos con y sin aisladores de base 



5 



Los porticos que contienen a la mamposteria son muy rigidos en comparacion con los 
porticos que no tienen mamposteria. Por lo tanto en cada direccion de analisis se tienen dos 
zonas, una fuerte y una debil; la primera se encuentra donde estan los porticos con la 
mamposteria acoplada y la segunda en el sentido opuesto que no tiene mamposteria. 

Para plantas muy alargadas, tanto los desplazamientos del portico fuerte como los del 
portico debil son mayores a los desplazamientos del Centro de Masas, CM. Ahora para 
estructuras torsionalmente rigidas, los desplazamientos laterales del portico debil son mayores 
a los desplazamientos del portico fuerte, en estructuras torsionalmente flexibles hay casos en 
que los desplazamientos del portico fuerte son mayores a los desplazamientos del portico debil. 
Hernandez y Lopez (2003). 

Lo cierto es que en un mismo sentido, los desplazamientos de los porticos exteriores 
son diferentes. Esta diferencia de desplazamientos es la que genera problemas de torsion en 
planta con resultados calamitosos que han llevado al colapso a las estructuras. Hernandez y 
Lopez (2007). 

Si no se considera la mamposteria de los porticos perimetrales y si la estructura es 
completamente simetrica, el Centro de Rigidez CR., coincide con el CM. En este caso los 
desplazamientos laterales de todos los porticos, en un sentido, son iguales. Al no coincidir el 
CR., con el CM., el desplazamiento lateral de los porticos es diferente y existen problemas de 
torsion. 



Se define el CR., como el lugar geometrico en el cual al aplicar las fuerzas sismicas la 
estructura se desplaza y no rota con respecto al piso inferior. El CR., solo existe en estructuras 
de un piso y en estructuras de varios pisos que sean compensables. Vasquez y Ridell (1984). 
En otro tipo de estructuras el CR., no existe en forma exacta pero existen procedimientos 
aproximados para determinar el CR., en forma aproximada. Cheung y Tso (1986), Villafane y 
Crisafulli (1986), Seguin (2007), entre otros. 

En este articulo se utiliza el procedimiento de Seguin (2007) para hallar el CR., en 
estructuras de varios pisos. Seguin (2007) de una estructura de multiples grados de libertad de 
varios pisos obtiene una estructura equivalente de un piso y con tres grados de libertad 
igualando las frecuencias y modos de vibracion. Es en esta estructura equivalente donde se 
halla el CR., se obtiene un solo CR., para todo el edificio. En Aguiar et al (2009) esta descrito 
el procedimiento de calculo del CR. 

Las estructuras construidas en la forma tradicional sin aisladores de base, en las que el 
CM., esta distante del CR., han colapsado o han tenido gran dano en sismos pasados. Ahora 
bien, las estructuras construidas sobre aisladores de base con mayores problemas de torsion 
se han comportado satisfactoriamente durante sismos recientes, como es el caso del Sur de 
California (no es un edificio abierto) pero si tiene una planta muy irregular en planta y se 
comporto en forma satisfactoria sin ningun dano en el sismo de Northridge de 1994, gracias a 
los aisladores de base elastomericos. Aguiar et al (2008). 



2. ANALISIS SISMICO EN ESTRUCTURAS SIN AISLADORES 



Dada una estructura de n pisos con base empotrada, sin aisladores, se encuentra las 
matrices de rigidez y de masas en un modelo de tres grados de libertad por planta. Con estas 
matrices se hallan los valores y vectores propios y posteriormente se encuentra el sistema 
equivalente de 1 piso con tres grados de libertad. Esquematicamente esto se ilustra en la 
figura 3, a la izquierda se muestra un edificio, que se analiza con tres grados de libertad por 
planta (losa rfgida) y a la derecha el modelo equivalente de un piso con tres grados de libertad, 
es en este modelo en que se hallan las excentricidades estaticas. 
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Ahora bien si se considera que el sismo actua en un solo sentido es factible analizarlo 
con dos grados de libertad como se indica en la figura 4. Donde e x es la excentricidad estatica 
en sentido X; u x (y) es la componente de desplazamiento horizontal en sentido Y; w,(0) es la 
rotacion de la planta con respecto al eje perpendicular de la losa. El sismo actua en sentido Y. 




Figura 4 Modelo numerico de calculo de una estructura monosimetrica. 



Notese que las dimensiones de la planta son 2a 1 2c. Siendo c<a. El sistema de 
ecuaciones diferenciales para el modelo de analisis de la figura 3, es: 
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Siendo m la masa equivalente del sistema de un piso, la misma que es igual a la 
suma de las masas de cada uno de los pisos de la estructura de n piso; K es la rigidez 
equivalente traslacional del sistema de un piso; p es el radio de giro del momento de inercia 
de la masa que se encuentra con la siguiente expresion. 



p = 2* 



2 . 2 

id + C 



12 



( 2 ) 



Antes de continuar indicando el significado de las restantes variables, es importante 
destacar que el un grado de libertad es el desplazamiento en sentido Y que se ha denominado 
u 1 (y)y el otro grado de libertad es p u x (0), de esta manera las dos coordenadas tienen las 
mismas unidades y son comprables. 



Las variables, no definidas todavia de ecuacion ( 1 ) son: e es la excentricidad estatica 
normalizada; Cl g es la relacion de la frecuencia torsional con respecto a la frecuencia 

traslacional; U es la aceleracion del suelo en sentido Y. 



P 



I V CR 

a e =:^— 

p* Kyy 



( 3 ) 

( 4 ) 



CR 

Donde K ee , es la rigidez a torsion pero referida al CR. La misma que se puede 
hallar a partir de K ee que es con respecto al CM. 



K 



CR 



ee 



= K 



ee 



K 



YY 



( 5 ) 



Siendo e x , la excentricidad estatica en sentido X. KyY,K ee son elementos de la 
diagonal de la matriz de rigidez del modelo de un piso y de dos grados de libertad. 

La matriz de amortiguamiento C (s) que se indica en la ecuacion ( 1 ) se halla con la 
siguiente ecuacion. 

C (s) = (6) 

Donde <1> es la matriz modal, conformada por los dos modos de vibracion, M es la 
matriz de masas indicada en la ecuacion ( 1 ). La matriz C es diagonal y vale: 
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C = 



2 £W X 



2%W 2 



(7) 



Siendo E, el factor de amortiguamiento de la estructura, normalmente se considera 
E, = 0.05 y es igual en todos los modos de vibracion. W 2 son las frecuencias de vibracion 
del modelo de 2 gdl. 

3. ANALISIS SISMICO EN ESTRUCTURAS CON AISLADORES 

Ahora, se considera que la estructura se halla sobre aisladores de base elastomericos, 
sin nucleo de plomo. Aguiar et al (2008). En este caso el modelo numerico de calculo es el 
indicado en la figura 5. Se ha notado con la letra u las coordenadas de la superestructura y 
con la letra q del sistema de aislacion. 




Figura 5 Modelo numerico de cuatro grados de libertad. 

Para el modelo de 4 grados de libertad, el sistema de ecuaciones diferenciales esta 
definido por: 
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Donde M (i) es la matriz de masas considerando la base empotrada del modelo de 2 
gdl; M <n es la matriz de masa total que se halla con la siguiente expresion: 

M (t) =r (s)t M (s) r (s) +M (b) (9) 

Siendo M u,) la matriz de masa del sistema de aislacion; C (s) , C ib) , son las matrices 
de amortiguamiento de la superestructura y del sistema de aislacion; K <s) , K ih) son las 
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matrices de rigidez de la superestructura y del sistema de aislacion. r is) es una matriz que se 
obtiene midiendo los grados de libertad en la superestructura correspondientes a las 

deformadas elementales del sistema de aislacion, para el presente caso r (s) = I . Por otro lado 
r ib) es el vector que se obtiene al medir los grados de libertad del sistema de aislacion si el 
suelo se desplaza la unidad, r a ” = [l 0] ' . Aguiar et al (2008). 



4. ESTRUCTURAS Y SISMOS DE ANALISIS 



Se analizaron estructuras de tres, cuatro, cinco y seis pisos. Para cada una de ellas se 
considero cinco plantas diferentes que son las indicadas en la figura 6, de tal manera que se 
tienen 20 edificios analizados sin aisladores de base e igual numero de edificios con aisladores 
de base. Cada uno de estos edificios fue sometido a la accion de 28 sismos artificiales. 




El parametro a relaciona la luz corta con respecto a la luz larga; este parametro varfa 
desde 1, que corresponde a una planta cuadrada hasta 0.5 que corresponde a una planta 
rectangular en la cual la luz corta es la mitad de la luz larga. La mampostena acoplada a los 
porticos se encuentra en el portico izquierdo (portico A) y en el portico superior (portico 1). Por 
lo tanto, son estos los porticos fuertes y los opuestos los porticos debiles. En la figura 6 se 
indica tambien la direccion en que actua el sismo. (Sentido Y). 

En Aguiar et al (2009) se indica las secciones de las columnas y vigas, las cargas 
actuantes; en todos los edificios la altura del primer piso es de 4.0 m., y la altura de los pisos 
superiores es de 3.0 m. Se considero que el hormigon utilizado tiene una resistencia a la 

compresion / c = 210 kg / cm 2 . La resistencia de la mampostena a la compresion 
f m =35 kg / cm 2 . El modulo de elasticidad del hormigon se obtuvo con E = 12000 ; el 



Roberto Aguiar Falcon!, Alex Jara y Dario Paez 



10 



modulo de elasticidad de la mampostena E m = 500 f m . El espesor de las paredes 
t = 0.15 m. 



En la figura 7 se presenta a la izquierda, los espectros de respuesta elasticos de cada 
uno de los siete sismos artificiales para suelo SI, se los ha identificado como 1S1, hasta 7S1. 
La duracion del 1S1 es 20 segundos y la duracion de los restantes sismos se incrementando 
en 5 s., en todos los sismos la fase intensa tiene una duracion de 10 s. En la figura izquierda 
tambien se ha dibujado con linea mas gruesa el espectro promedio. 



A la derecha de la figura 7 se muestra el espectro promedio, de los siete sismos 
artificiales y el espectro del CEC-2000 para el perfil de suelo SI . Se observa que los dos son 
parecidos ya que los sismos artificiales generados son compatibles con la forma espectral del 
CEC.2000. 




— 3S1 
— 4S1 




Figura 7 Espectros de acelerogramas artificiales, espectro promedio y del CEC-2000. Suelo SI . 



En el CEC-2000 se definen 4 tipos de suelo identificados por SI (suelo de resistencia 
muy dura); S2 (suelo de resistencia intermedia); S3 (suelo blando) y S4 (suelo muy blando). 
Para cada uno de estos perfiles de suelo se generaron 7 sismos artificiales, dando un total de 
28 sismos con los que se trabajo. 



5. RESULTADOS OBTENIDOS 



En esta apartado se presentan las caractertsticas del modelo equivalente de un piso 
con tres grados de libertad de cada uno de los edificios estudiados. Posterioremente, para uno 
de ellos se encuentra la respuesta en el tiempo con y sin aisladores de base con el proposito 
de que se visualice el trabajo realizado. Finalmente, se hallan los factores de amplificacion por 
torsion de los porticos exteriores, los mismos que fueron evaluados con las siguientes 
ecuaciones para las estructuras sin aisladores de base. 



r = U- a ~U CM 



r. = — 



U CM 

u - u 



CM 



U 



CM 



( 10 ) 



Donde r_ a ,r o son los factores en el portico fuerte (portico A) y en el portico debil 
(portico opuesto al A); U_ a ,U a son los desplazamientos laterales en los porticos extremos 
(fuerte y debil); U CM es el desplazamiento lateral en el CM., para a - l serfa en el portico B. 
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Para las estructuras con aisladores de base, los factores de amplificacion de la 
superestructura se encuentran con las siguientes expresiones. 



Mb) 



p(i>) _ 



u™ - u 



(b) 

CM 



u 



(b) 

—a 



u 



C b ) 

CM 



u w -u 



( b ) 
CM 



u 

u 



—a 

(b) 



(11 ) 



u 



0 b ) 
CM 



u„ 



Siendo , r^' 1 los factores de amplificacion por torsion en la superestructura con 
aisladores de base; , U ( a b) , U ( c b ^ los desplazamientos en la superestructura con aisladores 

de base en el portico fuerte, debil y CM; U _ a , U a los desplazamientos en los porticos fuerte y 
debil de la estructura con base empotrada, sin aisladores. 

De esta manera se incorpora en el calculo del factor de amplificacion por torsion la 
reduccion en desplazamientos de la estructura con aisladores de base en relacion a los 
desplazamientos con base empotrada. 



5.1 Rigidez y excentricidad de estructuras analizadas 



En la figura 8 muestra la ubicacion del CR., en las estructuras de 3 a 6 pisos para una 
planta rectangular con a - 0.5 . Como era de esperarse los CR. se halla a la izquierda ( e x ) y 
en la parte superior ( e Y ) del CM. Esto es para el caso de la estructura con base empotrada. 



Portico Fuerte 




Figura 8 Ubicacion del CR., en edificios abiertos con Ci — 0.5 
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En la tabla 1, se indican los valores de las excentricidades estaticas para cada uno de 
los edificios analizados, los elementos de la diagonal de la matriz de rigidez del modelo 
equivalente de 3 gdl., el momento de inercia de las masas y la frecuencia torsional. 



Tabla 1 Propiedades dinamicas de estructuras analizadas. 



3 PISOS 



Relacion de 
luz de losa 


ex (m.) 


ey (m.) 


Kxx 


Kyy 


K00 


P 


Q0x 


O0y 


a = 1 


-2.980 


2.980 


6150 


6150 


297890 


4.08 


2.110 


1.320 


a = 0.67 


-3.516 


2.774 


11150 


8710 


628330 


5.2 


1.730 


1.250 


a = 0.5 


-4.211 


2.843 


6710 


4710 


910390 


6.45 


1.660 


1.166 


a = 0.625 


-3.385 


2.428 


13635 


9876.6 


745540 


5.45 


1.664 


1.227 


a = 0.83 


-2.839 


2.363 


10367 


8746.2 


465350 


4.51 


1.719 


1.302 



4 PISOS 



Relacion de 
luz de losa 


ex (m.) 


ey (m.) 


Kxx 


Kyy 


K00 


P 


O0x 


O0y 


a = 1 


-2.598 


2.598 


5960 


5960 


270840 


4.08 


1.920 


1.350 


a = 0.67 


-2.820 


2.301 


10910 


8770 


595260 


5.2 


1.659 


1.302 


a = 0.5 


-3.329 


2.367 


14380 


10110 


989040 


6.45 


1.586 


1.210 


a = 0.625 


-2.668 


1.954 


13834 


10132 


720810 


5.45 


1.954 


1.258 


a = 0.83 


-2.275 


1.893 


10533 


8795.3 


448050 


4.51 


1.650 


1.343 



5 PISOS 



Relacion de 
luz de losa 


ex (m.) 


ey (m.) 


Kxx 


Kyy 


K00 


P 


O0x 


O0y 


a = 1 


-2.166 


2.166 


6020 


6020 


743860 


4.08 


1.803 


1.404 


a = 0.67 


-2.230 


1.873 


11140 


9180 


590580 


5.2 


1.596 


1.331 


a = 0.5 


-2.603 


1.939 


14670 


10760 


988860 


6.45 


1.524 


1.234 


a = 0.625 


-2.071 


1.554 


14566 


10782 


728440 


5.45 


1.540 


1.268 


a = 0.83 


-1.786 


1.493 


11116 


8746.2 


465350 


4.51 


1.595 


1.361 



6 PISOS 



Relacion de 
luz de losa 


ex (m.) 


ey (m.) 


Kxx 


Kyy 


K00 


P 


Q0x 


O0y 


a = 1 


-1.679 


1.679 


6150 


6150 


259800 


4.08 


1.763 


1.454 


a = 0.67 


-1.772 


1.518 


11760 


9900 


606550 


5.2 


1.542 


1.340 


a = 0.5 


-2.040 


1.578 


15500 


11800 


1023900 


6.45 


1.472 


1.243 


a = 0.625 


-1.617 


1.230 


15799 


12759 


760580 


5.45 


1.497 


1.263 


a = 0.83 


-1.406 


1.178 


12054 


9911.1 


469430 


4.51 


1.552 


1.360 



5.2 Respuesta en el tiempo 



Para ilustrar el trabajo realizado y para entender el problema de la torsion se presenta 
la respuesta en el tiempo de una estructura de 3 pisos con planta rectangular de dimensiones 
10 m., por 20 m., ( a = 0.5). Primero se presentan los resultados con base empotrada y luego 
con aisladores de base. 
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Para la estructura con base empotrada, los datos son los siguientes: e x = 4.211m 
excentricidad estatica; e x =0.653 excentricidad estatica normalizada (ecuacion 3 ); 
T s = 0.474 s perfodo de vibracion en sentido Y.; m = 36.73 T s 2 / m la masa equivalente 
que es igual a la suma de las masas de cada piso; = 1.166 que es la frecuencia torsional 
en sentido Y; E, ys) = 0.05 factor de amortiguamiento de la estructura. 

Se encuentra la respuesta en el tiempo, para el acelerograma 1S2 que tiene una 
duracion de 20 s., y corresponde a un perfil de suelo S2. Las respuestas se presentan en las 
figuras 9 y 10. 




Figura 9 Respuesta en el tiempo de estructura con base empotrada. 



En la parte superior izquierda de la figura 9 se muestra la historia de desplazamientos 
en el CM., de la estructura. Notese que estos desplazamientos se encuentran entre 5 y 8 cm. 
En la figura 10 se vera como se amplifican los desplazamientos en los porticos exteriores por 
efecto de la torsion. 

En la parte superior derecha de la figura 9 se muestra el histograma de la excentricidad 
dinamica. La excentricidad dinamica se obtiene dividiendo el momento de torsion con respecto 
al cortante. Se aprecia en el eje de las X los valores de esta excentricidad y se ve que no son 
cero por que existe torsion. Lo ideal es que en el histograma se tenga todos los valores 
concentrados en el punto cero pero esto se da en estructuras que no tienen problemas de 
torsion. 
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En la grafica inferior de la figura 9, se presenta la repuesta del giro de torsion 
normalizado y del desplazamiento lateral, cada punto representa como se desplaza y rota la 
estructura. 




Figura 10 Historia de desplazamientos en porticos exteriores sin aisladores de base. 



En la figura 10 se muestra la respuesta en desplazamientos de los porticos debil, en la 
parte superior y fuerte en la parte inferior. Se aprecia que estos desplazamientos son mayores 
a los desplazamientos en el CM., debido a la torsion. Los desplazamientos maximos del portico 
debil estan alrededor de los 30 cm., cantidad muy alta comparando con los desplazamientos 
del CM. En el portico fuerte estos desplazamientos maximos estan en por los 16 cm. 

Por otra parte, se han colocado aisladores de base de tal manera que el perfodo 
objetivo T ih) =2.0 s, que el CR coincida con el CM., luego e a,) =0. La relacion de 
frecuencia a torsion con respecto a la frecuencia traslacional =1.2649. El factor de 
amortiguamiento del sistema de aislacion £ ib) =0.10; la masa del sistema de aislacion se 

considera el 10% de la masa de la superestructura, m {b) = 3.673 T s 2 / m. La respuesta en 
desplazamientos del sistema de aislacion y de la superestructura se muestra en la figura 1 1 . 



Al coincidir el CM., con el CR., en el aislamiento, todos los porticos del sistema de 
aislamiento que estan en una misma direccion se desplazan lo mismo. En la superestructura 
continua habiendo torsion ya que los porticos extremos no son iguales como se aprecia en la 
parte inferior de la figura 11 sin embargo estos desplazamientos son mucho mas bajos a los 
que se tenia en la estructura con base empotrada. 
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Figura 12 Desplazamientos en porticos exteriores de superestructura con y sin aisladores. 
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En efecto, en la figura 12 se muestran los desplazamientos laterales de los porticos 
exteriores de la superestructura sin y con sistema de aislamiento, se aprecia que los 
desplazamientos de la estructura sin aisladores de base son mucho mas altos que los 
desplazamientos de la estructura con aisladores de base, con lo que se reduce notablemente el 
problema de torsion ya que al ser los desplazamientos bajos los momentos de torsion que se 
van a generar son bajos. 




Figura 13 Desplazamientos en sistema de aislamiento e histograma de excentricidad dinamica. 



En la figura 13 se aprecia a la izquierda los desplazamientos del sistema de aislacion y 
a la derecha el histograma de la excentricidad dinamica, se aprecia que todos los valores estan 
en un solo rango que coincide el punto 0. Esto significa que no existe torsion en el sistema de 
aislacion. Este histograma es de desear que tengan las estructuras con y sin aislamiento. 



Tabla 2 Valores medios de desplazamientos obtenidos, para Ct — 0.5 y en suelo S2. 



PISOS 


Desplazamientos Estructura sin 
sistema de aislamiento 


Desplazamientos Estructura con 
sistema de aislamiento 


SISMOS 


Portico Debil 


CM 


Portico Fuerte 


Portico Debil 


CM 


Portico Fuerte 


Promedio 


0.3190 


0.0828 


0.1641 


0.0600 


0.0383 


0.0166 



4 PISOS 


Desplazamientos Estructura sin 
sistema de aislamiento 


Desplazamientos Estructura con 
sistema de aislamiento 


SISMOS 


Portico Debil 


CM 


Portico Fuerte 


Portico Debil 


CM 


Portico Fuerte 


Promedio 


0.3416 


0.0871 


0.1749 


0.0669 


0.0458 


0.0246 



5 PISOS 


Desplazamientos Estructura sin 
sistema de aislamiento 


Desplazamientos Estructura con 
sistema de aislamiento 


SISMOS 


Portico Debil 


CM 


Portico Fuerte 


Portico Debil 


CM 


Portico Fuerte 


Promedio 


0.3305 


0.0903 


0.1736 


0.0716 


0.0523 


0.0329 



6 PISOS 


Desplazamientos Estructura sin 
sistema de aislamiento 


Desplazamientos Estructura con 
sistema de aislamiento 


SISMOS 


Portico Debil 


CM 


Portico Fuerte 


Portico Debil 


CM 


Portico Fuerte 


Promedio 


0.3622 


0.1001 


0.1962 


0.0757 


0.0583 


0.0410 
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5.3 Valores promedios de respuestas con y sin aisladores para a = 0.5 



Se procedio en forma similar a la indicada en el apartado anterior con todas las 
estructuras y con todos los sismos. Luego de ello se obtuvieron las respuestas maximas y de 
estas se obtuvo el promedio. Para ilustrar lo indicado, en la tabla 2 se muestran los 
desplazamientos laterales encontrados en los porticos exteriores de las estructuras con y sin 
aisladores de base, para las estructuras que tiene a = 0.5 . Los valores indicados es el 
promedio de las respuestas maximas encontradas con los 7 sismos artificiales para suelo S2. 

En la tabla 2 se aprecia que en la estructura sin aisladores de base, los 
desplazamientos laterales en los porticos exteriores son mayores a los desplazamientos en el 
CM., con lo que el problema de la torsion es muy grave. Ahora, en la estructura con aisladores 
de base se aprecia que los desplazamientos laterales del portico debil son mayores a los del 
CM., pero los desplazamientos laterales del portico fuerte son menores, mejorando 
notablemente el comportamiento de la estructura. 



Tabla 3 Factores de amplificacion por torsion en estructura sin y con aisladores. 



3 PISOS 


Factor de Amplificacion 
Torsional sin Aislamiento 


Factor de Amplificacion 
Torsional con Aislamiento 


SISMOS 


Portico Debil 


Portico Fuerte 


Portico Debil 


Portico Fuerte 


Promedio 


2.8513 


0.9809 


0.1065 


-0.0574 



4 PISOS 


Factor de Amplificacion 
Torsional sin Aislamiento 


Factor de Amplificacion 
Torsional con Aislamiento 


SISMOS 


Portico Debil 


Portico Fuerte 


Portico Debil 


Portico Fuerte 


Promedio 


2.9202 


1.0071 


0.0905 


-0.0651 



5 PISOS 


Factor de Amplificacion 
Torsional sin Aislamiento 


Factor de Amplificacion 
Torsional con Aislamiento 


SISMOS 


Portico Debil 


Portico Fuerte 


Portico Debil 


Portico Fuerte 


Promedio 


2.6578 


0.9220 


0.0797 


-0.0703 



6 PISOS 


Factor de Amplificacion 
Torsional sin Aislamiento 


Factor de Amplificacion 
Torsional con Aislamiento 


SISMOS 


Portico Debil 


Portico Fuerte 


Portico Debil 


Portico Fuerte 


Promedio 


2.6194 


0.9605 


0.0622 


-0.0621 



En la tabla 3 se presentan los factores de amplificacion por torsion, en los porticos 
exteriores para la estructura calculada sin aisladores de base, a la izquierda y con aisladores 
de base a la derecha. Cuando la estructura no tiene aisladores de base los factores de 
amplificacion son altos, en el portico debil vartan entre 2.6 y 2.9. Esto significa que las fuerzas 
laterales en los porticos debiles seran aproximadamente 3 veces mayores a las fuerzas 
laterales con las cuales fueron disenadas si no se considero la mamposteria. 



En cambio cuando se tiene aisladores de base estos factores de amplificacion por 
torsion, son bastante bajos con lo que practicamente se ha eliminado el problema de la torsion 
en edificios abiertos. 
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Figura 14 Factores de amplificacion en cada uno de los edificios analizados. 
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5.4 Factores de amplificacion en estructuras sin aisladores de base 



En la figura 14 se muestran los factores de amplificacion por torsion en las estructuras 
con base empotrada, para todos los casos analizados, para un perfil de suelo SI. La primera 
grafica corresponde a plantas rectangulares con a - 0.5 y la ultima que aparece en la parte 
inferior derecha para a - 1.0. En cada una de ellas se tienen los factores de amplificacion de 
desplazamientos para los edificios de 3 a 6 pisos. 



Mientras mas alargada es la planta, mayores son los factores de amplificacion, para el 
portico debil. 



6. FACTORES DE AMPLIFICACION 



En la figura 15 se muestran los factores de amplificacion encontrados en las 
estructuras con plantas cuadradas, de 3 a 6 pisos; en los diferentes tipos de suelo, con y sin 
aisladores de base. Se aprecia que el tipo de suelo no influye en el valor de T . Es verdad que 
se tienen diferentes desplazamientos por el tipo de suelo pero al calcular T se tienen valores 
similares en los cuatro tipos de suelo. 




Figura 15 Factores de amplificacion con y sin aisladores. Para CC = 1 



En la figura 16 se presentan los resultados para a = 0.83; en la figura 17 para 
a = 0.67 ; en la figura 1 8 para a = 0.625 y en la figura 1 9 para a = 0.5 . 

En todos los casos analizados se observo que el suelo no influye mayormente en el 
calculo del factor de amplificacion. A excepcion de la planta rectangular con a = 0.5 en que si 
se notan diferencias en el factor de amplificacion, especialmente para el portico debil. 



Recibido: Julio de 2008 
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Figura 16 Factores de amplificacion con y sin aisladores. Para CC = 0.83 





Figura 17 Factores de amplificacion con y sin aisladores. Para CC — 0.67 



Lo que si se nota en las figuras 15 a 19 es que se reduce notablemente el factor de 
amplificacion por torsion en las estructuras con aisladores de base. En la mayor parte de los 
casos se aprecia que este factor esta alrededor de cero, lo que significa que no existen 
problemas de torsion en las estructuras con aisladores de base. 



La reduccion de la torsion en estructuras con aisladores de base en gran parte se debe 
a que los desplazamientos laterales en la superestructura son muy bajos en comparacion con 
los que se producen en el sistema de aislamiento. Esto influye en que los valores de Y sean 
bajos al utilizar la ecuacion ( 11 ). 
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7. CONCLUSIONES 



Se ha encontrado el factor de amplificacion por torsion en 20 edificios abiertos de 
hormigon armado de 3 a 6 pisos, con diferentes dimensiones en planta, ante la accion de 28 
sismos sinteticos compatibles con los espectros del Codigo Ecuatoriano de la Construccion, 
para la zona de mayor peligrosidad stsmica del Ecuador y en los cuatro perfiles del suelo que 
contempla el codigo. Se han analizado los edificios con y sin aisladores de base. Del estudio 
realizado se despenden las siguientes conclusiones: 

• En los edificios abiertos con base empotrada, la excentricidad estatica esta alrededor 
del 20 al 30% de la longitud de la planta, lo que implica que tienen importantes 
problemas de torsion. 
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• El problema de torsion es mas grave en edificios con planta rectangular en que la una 
dimension es el doble de la otra, a - 0.5 . En estos edificios el factor de amplificacion 
por torsion esta alrededor de 3 para el portico debil y en el portico fuerte los 
desplazamientos laterales son mayores a los desplazamientos del CM. De tal manera 
que es muy critico el problema en este tipo de estructuras con a = 0.5 . En edificios 
abiertos con planta cuadrada tambien existen problemas de torsion pero el problema 
es menos critico. 

• Los aisladores de base minimizan los problemas de torsion en edificios abiertos, en 
parte debido a que los desplazamientos laterales en la superestructura son bastante 
bajos en comparacion con los desplazamientos que se producen en el sistema de 
aislamiento. Esto se da cuando en el sistema de aislacion se tiene que el centro de 
masas coincide con el centro de rigidez es decir no hay excentricidad estatica en el 
sistema de aislamiento. Por lo tanto, si se desea construir un edificio abierto en el cual 
la mamposteria se acopla a la estructura lo ideal es que este se construya sobre 
aisladores de base. 
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RESUMEN 

La amplificacion directa de DNA bacteriano a partir de muestras clinicas con cebadores 
de PCR de amplio espectro es una importante alternativa para el reconocimiento de patogenos 
en sangre. Sin embargo, la principal desventaja de esta tecnica es la dificultad para preparar la 
muestra clfnica para el ensayo. 

El proposito de este estudio fue disenar un metodo de alta capacidad de procesamiento 
de muestras con eliminacion del background genomico humano. Con este objetivo, se modified 
un sistema de extraccion de DNA basado en el uso de columnas con membranas de silica gel, 
usando un paso previo de centrifugacion diferencial. La sensibilidad analitica de este metodo 
de purificacion de DNA fue evaluada con un ensayo 16S rDNA PCR mediante el uso de 
diluciones seriadas de cultivos bacterianos de especies modelo. 

El background genomico humano fue eliminado por remocion del plasma y capa leuco- 
plaquetaria. Se comprobo que este paso no reduce la cantidad de CFU presentes en la 
muestra. El metodo disenado demostro ser eficiente y reproducible. La sensibilidad del metodo 
fue 1 CFU/mL. Ademas, probo ser aplicable para el uso con muestras clinicas. Este estudio 
sera complementado con una investigacion donde se validara clfnicamente el ensayo 
descontaminado con muestras clinicas reales de pacientes con bacteriemia/sepsis. 

ABSTRACT 

Direct amplification of bacterial DNA from clinical specimens with broad-range primers 
through PCR provides an important alternative to the detection and identification of pathogens 
causing bacteremia/sepsis. Nevertheless, the main drawback of this method is the difficulties in 
sample processing to comply with the requirements of DNA extraction systems in terms of 
sample volume and DNA extraction capacity. 

The purpose of this study was to design a high capacity sample processing method with 
elimination of background human DNA while conserving pathogen-related DNA. For this goal, a 
DNA extraction procedure based in a silica gel membrane column was carried out including a 
previous differential centrifugation step. The analytical sensitivity of this DNA purification method 
was evaluated by a 16S rDNA PCR on dilution series of model bacteria spiked into blood 
samples. 

The background human DNA was eliminated by removal of plasma and buffy coat. It 
was demonstrated that this step does not decrease the amount (CFU) of bacteria present in the 
sample. The method optimized here showed to be efficient and reproducible. With the previous 
removal of human DNA from the sample, the sensitivity of the PCR assay reached 1 CFU/mL. 
We applied the assay for the detection of pathogenic DNA in 2 samples of a hemoculture- 
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positive bacteremic patient with excellent results. The present study will be complemented by a 
clinical validation step by applying our decontaminated PCR assay on clinical samples from 
patients with clinically suspected bacteremia/sepsis. 

1 INTRODUCCION 

La infeccion y sepsis son problemas clinicos especialmente importantes, en donde la 
aplicacion de nuevas tecnicas diagnosticas podrfa contribuir en terminos de disminucion de 
tasas de morbi-mortalidad, reduccion de estancia hospitalaria, decremento de incidencia de 
infeccion nosocomial y disminucion global de los costos clinicos. 

La mayor parte de las bacterias de interes clfnico pueden identificarse mediante 
tecnicas microbiologicas convencionales, las cuales precisan del aislamiento previo del 
patogeno y se fundamentan en sus propiedades fenotfpicas (Rocidio y Mendoza, 2004). 

El hemocultivo se ha convertido en el estandar de oro para el diagnostico de sepsis. 
Sin embargo, esta tecnica puede no ser lo suficientemente rapida debido a que usualmente 
toma entre 2 - 8 dfas para proporcionar resultados definitivos (Crowley, 2002). 

En los ultimos anos se ha promovido el uso de tecnicas moleculares como PCR 
(Polymerase Chain Reaction) y secuenciacion para el diagnostico de enfermedades causadas 
por bacterias. La amplificacion directa del DNA bacteriano a partir de muestras clfnicas 
mediante el uso de PCR de amplio espectro presenta ventajas para la determinacion de 
patogenos infecciosos en tejidos comunmente esteriles (Rantakokko-Jalava & Jalava, 2002). 

El aislamiento de DNA es una tecnica basica en biologfa molecular. La calidad, 
cantidad e integridad del DNA afectan directamente a los resultados de las pruebas efectuadas 
a partir de este. Segun Surzycki (2000), un procedimiento de extraccion a gran escala puede 
proporcionar de 200 - 300 pg de DNA cromosomal a partir de 2 - 4 x 1 0 10 bacterias frescas o de 
20 - 30 mg de biomasa congelada. No obstante, cuando el DNA bacteriano proviene de 
muestras clfnicas, su probabilidad de recuperacion se reduce conforme disminuye el volumen 
de muestra. Al mismo tiempo, el background genomico humano se convierte en una barrera 
ffsica para la adhesion del DNA patogeno a diversas membranas o partfculas empleadas para 
la purificacion, haciendo el ensayo poco reproducible (Jordan & Durso, 2000). 

La variabilidad de tipos de muestras clfnicas incrementa la complejidad del diagnostico 
bacteriano por PCR. Esto se ve reflejado en la dificultad de extrapolar los Ifmites de deteccion 
determinados experimentalmente a la verdadera sensibilidad en muestras clfnicas 
(Rantakokko-Jalava & Jalava, 2002). 

Muchos investigadores han planteado la posibilidad de que el background genomico 
humano sea un factor inhibitorio durante reacciones PCR (Rantakokko-Jalava et at, 2002). Un 
metodo optimo de procesamiento de muestras, por lo tanto, debera concentrar el DNA, 
fundamentalmente el perteneciente al organismo bianco y eliminar los agentes inhibitorios 
usualmente presentes en los fluidos biologicos. 

A pesar de la amplia utilizacion de ensayos moleculares para estudios clinicos, no 
existen investigaciones a gran escala de validacion de sistemas moleculares en este escenario 
especffico. Un taxativo importante ha sido la falta de estandarizacion y reproducibilidad en el 
uso de la tecnica de extraccion de DNA de patogenos en muestras clfnicas, asf como la 
carencia de sistemas de extraccion de alta capacidad que soporten grandes volumenes de 
muestra. 

Este trabajo propone la modificacion de un sistema de extraccion y purificacion de DNA 
de patogenos presentes en muestras clfnicas mediante la eliminacion de DNA genomico 
humano de leucocito. Mediante esta modificacion se espera obtener un sistema capaz de 
tolerar altos volumenes de muestra sin que la capacidad de extraccion se vea alterada. 
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2 MATERIALES Y METODOS 

2.1 Cultivos puros bacterianos 

Se realizaron siembras microbiologicas en agar a partir de cultivos bacterianos puros 
de Escherichia coli, Staphylococcus aureus y Enterococcus faecaiis, (provistos por I. Narvaez, 
Hospital Carlos Andrade Marin), con el fin de obtener el volumen de bacterias necesario para la 
transferencia a medio de cultivo liquido. Los medios de cultivo empleados fueron escogidos en 
base a los requerimientos de crecimiento de cada bacteria. Se utilizaron los medios LEVINE 
EMB Agar ( Escherichia coli) (ACUMEDIA, U.S.A.), Manitol Salt Phenol-Red Agar (ACUMEDIA, 
U.S.A.) (Staphylococcus aureus ) y Agar Sangre (BECTON DICKINSON, Francia) 
(Enterococcus faecaiis). 

Posteriormente se inocularon tubos con LB Broth (BECTON DICKINSON, Francia) para 
conseguir el crecimiento exponencial de cada tipo de bacteria. 

En cada siembra se colocaron tres asas colmadas de bacterias (aproximadamente 150 
pL), tomadas de los cultivos bacterianos en agar. El experimento se hizo por duplicado. Los 
tubos inoculados fueron llevados a incubacion con agitacion a 37 S C durante 18h. 

2.2 Preparacion de Diluciones Seriadas 

Se realizaron diluciones seriadas en suero fisiologico a partir de los cultivos puros 

obtenidos en el paso descrito previamente. Las diluciones obtenidas fueron valoradas mediante 
un ensayo de conteo de unidades formadoras de colonias (CFU). Los medios de cultivo 
utilizados para cada cepa bacteriana fueron los detallados anteriormente. 

2.3 Preparacion de Muestras Artificiales 

Para desarrollar un sistema de extraccion de DNA de patogenos a partir de sangre 
periferica fue necesario construir muestras infectadas con una cantidad de bacterias conocida 
para el posterior analisis de la tecnica. Para ello se anadio un volumen determinado de dilucion 
seriada de cultivo bacteriano en un volumen establecido de sangre completa (1:10). Las 
muestras artificiales fueron evaluadas mediante un ensayo de conteo de CFU. 

2.4 Eliminacion de DNA Genomico Humano 

Se probo la tecnica de centrifugacion diferencial (CD) con el objetivo de eliminar el DNA 
humano presente en leucocitos mediante la remocion del plasma y la capa leuco-plaquetaria. 
Las muestras fue sometida a este procedimiento antes de la extraccion de DNA y los 
resultados fueron evaluados en relacion a muestras no tratadas. 

Las muestras escogidas para CD fueron sometidas a centrifugacion por 10 minutos a 
3000 rpm, seguidas por eliminacion del plasma y la capa leuco-plaquetaria. Las muestras asi 
tratadas fueron posteriormente evaluadas mediante conteo de CFU. 



2.5 Extraccion de DNA a partir de Muestras de Sangre 

Se probaron dos metodos de extraccion de DNA que probaron ser eficientes para el 
aislamiento de acidos nucleicos de bacterias Gram-positivas y Gram-negativas. El primer 
protocolo fue una modificacion del kit QIAGEN DNA Blood Mini Kit (QIAGEN, Alemania) para 
adaptarlo para extraccion de DNA genomico bacteriano. El segundo procedimiento se baso en 
extraccion organica mediante cloroformo-isoamil alcohol La eleccion del mejor sistema de 
purificacion de DNA dependio de su eficiencia para la extraccion simultanea de acido nucleico 
bacteriano y humano. 
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Los sistemas empleados se basan en el uso de un coctel litico conformado por 
lisozima, lisostafina y liticasa, a concentraciones de 15mg/ml, 1.8mg/ml y 5U/pl 
respectivamente, para la destruccion de la pared celular bacteriana. La verificacion de la 
existencia de DNA en las muestras se realizo mediante electroforesis en gel de agarosa al 
0 . 8 %. 



2.6 Identificacion de DNA Genomico Bacteriano mediante 16S rDNA PCR 

Para la identificacion de la presencia de DNA genomico bacteriano se decidio emplear 
un ensayo de PCR de amplio espectro que amplifica un fragmento del gen rrs (16S rDNA). Las 
secuencias de los cebadores empleadas fueron tomadas de Warwick, et al (2003) y el 
protocolo de PCR de Izquierdo (2007). Los cebadores Nf (5’ TCC TAG GGG AGG CAG CAG T 
3’) y Nr (5’ GGA CTA CCA GGG TAT CTA ATC CTG TT 3’) poseen una longitud de 1 9 y 26 bp, 
respectivamente, y amplifican un fragmento de 476 bp. 

La visualizacion del resultado de la amplificacion por PCR se realizo mediante 
electroforesis en gel de agarosa al 2.0%. Con esta tecnica se verified la existencia del amplicon 
deseado. 

3 RESULTADOS Y DISCUSION 

Debido a que los agentes etiologicos mas comunes de bacteriemias en neonatos son 
Escherichia coli (Gram-negativa) y Staphylococcus aureus (Gram-positiva), se decidio 
emplearlos como microorganismos modelo de este estudio (Reimer et al, 1997). Sin embargo, 
S. aureus no se comporto adecuadamente. Se presume que este hecho se debe a que, las 
cepas empleadas para los ensayos fueron sometidas a un estres excesivo y por lo tanto 
perdieron su viabilidad. Para superar este inconveniente, se decidio emplear otra cepa 
bacteriana Gram-positiva, Enterococcus faecalis, con la que se obtuvieron los resultados 
deseados que se muestran en la figura 1. 

Conociendo de antemano que las bacterias pueden adherirse a las celulas sanguineas, 
en especial leucocitos, mediante fuerzas ffsico-quimicas o por especificidad del 
microorganismo a determinada superficie (Henderson et al, 2000), fue necesario el analisis de 
muestras contaminadas a traves de recuento de CFU, para de esta manera conocer si la 
perdida de bacterias en el proceso de eliminacion de componentes de la sangre afectaba de 
manera significativa los resultados del ensayo. 



Escherichia coli 



Enterococcus faecalis 





Figura 1 Resultados del recuento de CFU de E. coliy E. faecalis de los ensayos de siembra de 

diluciones seriadas. 



Para conseguir la eliminacion del DNA genomico humano, se decidio emplear la 
tecnica de centrifugacion diferencial, ya que con esta se consigue separar muestras 
sanguineas en sus componentes principales: plasma, leucocitos y eritrocitos (Anaya et al, 
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2005). La fase de interes para el estudio fue la capa eritrocitaria, por lo que al descartar el 
plasma y la capa leuco-plaquetaria se cumplio con el fin de anular el background genomico 
humano presente en linfocitos y mantener a las bacterias en la fraccion restante de la muestra. 
La evaluacion de este proceso demostro que no existe una perdida representativa de 
microorganismos, es decir que la eliminacion del DNA humano no afecta a los resultados que 
se puedan obtener a partir de la muestra. Ademas, se consigue la concentracion de las 
bacterias en una pequena fraccion de la sangre, hecho que aporta al sistema de aislamiento de 
DNA como se observa en la figura 2. 




Figura.2 Fotografia del recuento de CFU de las muestras contaminadas con E. faecalis tratadas con 
centrifugacion diferencial (CD): a) MEF8, b) MEF7 y c) MEF6 y sus controles: d) MEF8, e) MEF7 y f)MEF6 



Una vez cumplido el analisis de muestras artificiales, se procedio a probar metodos de 
aislamiento de DNA que provean el material genetico necesario para la evaluacion del sistema, 
que en este caso se refiere a la totalidad del DNA presente en la muestra artificialmente 
contaminada, ya sea humano o bacteriano. 

El posible organismo bianco de una prueba de PCR de amplio espectro puede ser 
altamente inestable en su capacidad de soportar procedimientos ffsicos y qufmicos 
(Rantakokko-Jalava & Jalava, 2002). Por lo tanto, se requiere que el metodo de aislamiento de 
DNA genomico bacteriano empleado pueda lisar las paredes mas resistentes (bacterias Gram 
positivas) sin danar indirectamente el DNA bianco purificado de especies mas fragiles 
(bacterias Gram negativas) (Coyne et al, 2004). Los dos sistemas examinados cumplieron con 
estas caracteristicas y estan basados en un coctel enzimatico con el fin de obtener DNA de 
ambos tipos de bacterias con igual eficiencia y sin comprometer a los acidos nucleicos de las 
especies mas sensibles a ruptura de la pared celular. 

El sistema de extraccion manual modificado no se acoplo a los requerimientos 
necesarios para la evaluacion del sistema debido a que se comporto selectivamente, es decir, 
no todas las especies de DNA en la muestra fueron aisladas con la misma eficiencia. Se 
presume que este hecho se debio a que la muestra no estuvo incubada por suficiente tiempo 
con proteinasa K, limitando la liberacion del DNA genomico humano. Sin embargo se comprobo 
con este metodo que la eliminacion del plasma y el buffy coat mejora la extraccion de DNA de 
patogenos en muestras clfnicas. 

El metodo modificado para aislamiento de acidos nucleicos por kit de adsorcion a 
columnas de silica fue altamente eficiente para la extraccion del DNA presente en la muestra 
(Figura 3). Por lo tanto, este sistema fue elegido para realizar la eliminacion de DNA genomico 
humano y posteriormente demostrar que la exclusion de este disminuye la interferencia que el 
DNA genomico produce en la reaccion de PCR y hace al ensayo mas reproducible. Ademas, 
este metodo fue usado para el procesamiento de altos volumenes de muestra sin perder su 
eficiencia. 
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1Kb M4 M5 M6 M4 M5 M6 CD4 CD5 CD6 CD4 CD5 CD6 CN 1Kb 



Figura 3 Electroforesis en gel de agarosa al 0.8% tenido con bromuro de etidio donde se observan los 
resultados del aislamiento de DNA con kit del ensayo preliminar del sistema disenado con empleo de un 
alto volumen de muestra. M4, M5 y M6 corresponden a las muestras artificiales de 1 mL sin tratar; CD4, 
CD5 y CD6 a las muestras de 5 mL tratadas por separacion de componentes de sangre y CN es un 

control negativo. 



El sistema de aislamiento disenado fue validado mediante el uso de 16S rDNA PCR 
(Figura 4). La amplificacion directa del DNA bacteriano a partir de muestras artificialmente 
contaminadas mediante el uso de PCR de amplio espectro es una ventaja para la 
determinacion de patogenos infecciosos en tejidos comunmente esteriles (Rantakokko-Jalava & 
Jalava, 2002). 




Figura 4 Electroforesis en gel de agarosa al 2% tenido con bromuro de etidio donde se observan los 
productos de PCR del ensayo preliminar del sistema disenado con empleo de un alto volumen de 
muestra. Los primeros dos carriles corresponden al marcador de peso molecular Low DNA Mass Ladder 
(INVITROGEN, U.S.A.), M4, M5 y M6 corresponden al DNA de las muestras artificiales de 1 mL sin tratar; 
CD4, CD5 y CD6 al DNA de las muestras de 5 mL tratadas por separacion de componentes de sangre y 

CN al control negativo. 



La presencia de un amplicon de 476 bp en el control negativo, en todos los ensayos de 
16S rDNA PCR, se genero por contaminacion de reactivos que provienen de una fuente 
bacteriana, como la Taq polimerasa. Este como otros reactivos puede contener DNA 
contaminante como resultado de su manufactura e incompleta purificacion (Corless et al 2000). 

El efecto de la eliminacion del background genomico humano en el momento del 
aislamiento se ve evidenciado en la Figuras 4. Se puede observar que los productos de PCR 
pertenecientes a las muestras procesadas sin ningun tratamiento adicional no presentan 
uniformidad alguna en la amplificacion. Mientras tanto, en las muestras tratadas por eliminacion 
del plasma y leucocitos se visualiza homogeneidad en los duplicados y entre muestras del 
mismo tipo. Estos resultados reflejan que el sistema disenado es reproducible y puede ser 
usado con grandes volumenes de muestra. 
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Con este sistema se puede detectar 1 CFU/mL de sangre (Figura 5). Por lo tanto, la 
implementacion de este metodo de aislamiento puede ayudar al diagnostico de bacteriemias de 
bajo nivel. 




Figura 5 Electroforesis en gel de agarosa al 2% tenido con bromuro de etidio donde se observan los 
productos de PCR del ensayo de determinacion del limite de deteccion del sistema de alta capacidad 
disenado para aislamiento de DNA de patogenos en muestras de sangre. El primer carril corresponde al 
marcadorde peso molecular Low DNA Mass Ladder (INVITROGEN, U.S.A, MC5, MC6, MC7 y MC8 
corresponden al DNA de las muestras artificiales de 1 mL sin tratar; CD5, CD6, CD7 y CD8 al DNA de las 
muestras de 5 mL tratadas por separacion de componentes de sangre y CN al control negativo. 



La evaluacion del sistema con muestras clfnicas, lo cual constituyo la ultima fase del 
presente estudio, fue exitosa. Se ratified el diagnostico por hemocultivo de sepsis y bacteriemia 
a partir de muestras de sangre periferica de una paciente del Hospital de Ninos Baca Ortiz de la 
ciudad de Quito mediante el uso de esta tecnica (Figura 6). 



1 (TObp- 1 1 CN 2 2 CN 




Figura 6 Electroforesis en gel de agarosa al 2% tenido con bromuro de etidio donde se observan los 
productos de PCR del ensayo de evaluacion del sistema de alta capacidad de aislamiento de DNA 
patogeno en sangre con muestras clfnicas. 1 y 2 corresponden a muestras clfnicas y CN a los controles 

negativos. 



Este estudio aporto a la superacion de los limitantes de las tecnicas de aislamiento de 
acidos nucleicos tradicionales mediante el diseno de una tecnica de extraccion de DNA que 
elimina el background genomico humano, incrementando la eficiencia de deteccion de 
patogenos en muestras clfnicas. Ademas, este metodo permite el procesamiento de altos 
volumenes de muestra, hecho que incrementa la sensibilidad del sistema lo cual podrfa facilitar 
el estudio de bacteriemias de bajo nivel. 
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RESUMEN 

El objetivo principal de este estudio fue caracterizar las fincas ganaderas en la zona de 
Santo Domingo de los Colorados y relacionar las variables socioeconomicas con su cobertura 
arborea. Para esto se realizaron encuestas al azar a 24 productores ganaderos de la zona, en 
donde se recopilo informacion acerca de las caracterfsticas de la finca, sistema de produccion 
ganadera, del componente arboreo y del componente herbaceo. El area total de las fincas 
encuestadas fue 2183 ha, de las cuales la mayor parte del area corresponds a potreros. Del 
total de fincas, el 29.17 % de las fincas fueron doble proposito, el 25 % mixtas (agricultura - 
ganaderfa), el 20 % de leche, el 16.67% de carne y el 8.33 agricolas. El 91 .67 % de las fincas 
presento arboles en potreros, la mayor cantidad de arboles en potreros reportados en la zona 
fueron usados como frutales y sombra para el ganado; las especies arboreas mas comunes en 
potreros mencionadas por los finqueros fueron Psidium guajaba, Citrus reticulata y Citrus 
sinensis. El 100 % de las fincas tuvieron cercas vivas, compuestas principalmente de Caraca 
(. Erythrina poeppigiana ); las especies presentes en cercas vivas fueron utilizadas como sombra 
y forraje para el ganado y como material vegetativo para propagacion. La cobertura arborea 
total (bosque primario + bosque secundario + area de plantaciones forestales + area de 
bosques riparios) de las fincas ganaderas estuvo influenciada en forma positiva por el 
porcentaje de area de cultivos y en forma negativa por el area de potreros, la produccion de 
leche, la edad del destete, el tiempo de pastoreo y los anos que los productores se han 
dedicados a la ganaderfa. 

Palabras claves: Fincas ganaderas; Santo Domingo de los Colorados; cobertura arborea; 
cercas vivas; arboles en potreros. 



1 INTRODUCCION 

En el Ecuador desde los anos cincuenta se viene observando una marcada 
disminucion de la cobertura forestal. Entre la decada de los setenta y hacia finales de los 
ochenta se noto un incremento en la deforestacion, llegandose a estimar tasas anuales de 
superiores a las 50,000 ha, debido al aumento de la frontera agrfcola, polfticas 
gubernamentales de titulacion de tierra, explotacion ilegal, incendios forestales, actividad 
ganadera extensiva, polfticas crediticias y desarrollo urbanfstico (Carrera, 1 999; Encarta, 2000). 

Muchas areas boscosas y las destinadas a la agricultura han sido modificadas a areas 
de pasturas debido a la expansion de la ganaderfa (Howard-Borjas 1995). Esta expansion ha 
estado asociada con la perdida de la sostenibilidad de los ecosistemas, la cual ha conllevado a 
la perdida y fragmentacion de bosques y la creacion de paisajes que son mosaicos de potreros, 
bosques y cultivos. Actualmente, en la Provincia de Pichincha existe una superficie de 421454 
hectareas destinadas a la actividad ganadera que corresponde al 45.5 % de la superficie total 
(MAG 2004). 

Ante esta situacion, el sector ganadero debera buscar alternativas de produccion que 
permitan producir competitiva y sosteniblemente, identificando en las fincas ganaderas las 
areas con mayor potencial productivo para cada alternativa (Ibrahim y Schlonvoigt 1999). Por 
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ello es importante conocer el estado actual del uso del suelo de las tineas ganaderas (bosques, 
cercas vivas, arboles en potreros, cultivos, plantaciones forestales), que puede ayudar al 
planeamiento y desarrollo de las estrategias a seguir para mejorar los estados de suelo 
existentes en beneficio de los sistemas productivos y de la conservacion del ecosistema; asi 
como tambien, es importante conocer los tipos de coberturas arboreas que se encuentran 
formando parte de los sistemas productivos y como estos varian entre las diferentes tineas en 
estructura, superficie que ocupan y composicion de especies; ademas de su utilidad, manejo e 
importancia como herramienta para la conservacion y la produccion (Harvey y Haber 1999). En 
este estudio se caracterizo la cobertura arborea en los diferentes tipos de sistemas productivos 
y se evaluo la relacion de la cobertura arborea con las variables socioeconomicas de los 
sistemas productivos. 



2 METODOS 

2.1 Localizacion y descripcion del area de estudio 



El presente trabajo se realizo en tineas ganaderas de Santo Domingo de los Colorados, 
Ecuador, area de influencia del proyecto CONTRIBUCION DEL COMPONENTE ARBOREO A 
LA PRODUCTIVIDAD Y CONSERVACION DE LA BIODIVERSIDAD EN FINCAS 
GANADERAS DE SANTO DOMINGO DE LOS COLORADOS. El area de estudio se encuentra 
en la provincia de Pichincha, Canton Santo Domingo de los Colorados, Region Litoral. La zona 
se ubica a una altura de 300 msnm, con una precipitacion media anual de 3.045.1 mm. La 
humedad relativa promedio es del 85.8 % y la temperatura promedio es de 24.2 S C con poca 
variacion durante el dia (INIAP 2004). 



2.2 Seleccion de las tineas 



Se realizo una consulta con productores de la zona, con la finalidad de identificar los 
diferentes estilos de vida de los productores y conocer las caracteristicas intrinsecas del area 
de estudio. Posteriormente se seleccionaron 24 tineas al azar de una lista de productores, 
proporcionada por la ASOGAN. Las tineas seleccionadas debfan cumplir con el requisito de 
estar activas. 

La recoleccion de datos se llevo a cabo por medio de encuestas estructuradas, con las 
cuales se recopilo informacion sobre las caracteristicas de cada tinea: uso del suelo (tamano y 
distribucion de las areas de la tinea); caracteristicas de los productores (edad, nivel de 
educacion; composicion y distribucion de la mano de obra); caracteristicas de los sistemas 
ganaderos (numero de cabezas, razas, suplementacion y tipos de pastos; caracteristicas, 
tamano y utilizacion de las areas agricolas y descripcion y composicion del componente 
arboreo (usos y especies de arboles en potreros y cercas vivas). 

2.3 Analisis de la informacion 

La informacion proveniente de encuestas de las 24 tineas fue analizada por medio de 
estadistica descriptiva (media, desviacion estandar, error estandar y coeficiente de variacion), y 
diferentes tecnicas de graficos. Para analizar las relaciones entre la cobertura arborea y las 
caracteristicas de las tineas se realizaron regresiones lineales y multiples entre las variables 
(INFOSTAT, 2002), de tal forma de identificar las variables socioeconomicas que tengan mas 
relacion con la cobertura arborea presente en las tineas. 
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3 RESULTADOS 



3.1 Caracterizacion general de las fincas en Santo Domingo de los Colorados, 
Pichincha en base a 24 encuestas 



3.1.1 Uso del suelo 

El area total de las 24 fincas encuestadas es 2183 ha, de las cuales la mayor parte del 
area correspond^ a potreros. En efecto, el 76 % de la zona fueron potreros, 15.4 % bosques 
primarios, 11.7 % a cultivos, 6.99 % bosques secundarios, 1.0 % pastos de code, 0.5 % 
plantaciones forestales, el resto corresponds a areas de construccion. 



3.1.2 Sistemas de produccion 

Los sistemas de produccion ganaderos en la zona de Santo Domingo de los Colorados 
son de diferentes tipos. El 29.17 % de las fincas son doble proposito, el 25 % se dedican a la 
combinacion agricultura - ganaderia, el 20 % son de leche, el 16.67% son de carne y el 8.33 
son agricolas. El tamano promedio de las fincas encuestadas fue de 90.96 + 17.29 ha, con 
variaciones entre 5 y 300 ha, lo que indica que hay mucha variabilidad entre el tamano de las 
fincas. 



La produccion bovina en el area de estudio vario notablemente entre las fincas. La 
carga animal fue de 10.63 UA/ha, con variaciones entre 0 a 180 UA/ha. El 87.5 % de las fincas 
ganaderas utiliza sistema de pastoreo rotacional con un perfodo de ocupacion entre 1 a 15 
dias y un periodo de descanso entre 28 y 90 dias. El resto de las fincas utilizan sistemas de 
pastoreo continuo. El tamano promedio de los potreros fue de 68.93 ± 12.94 ha, con un rango 
de 0.1 a 200 ha. Otros parametros productivos de las fincas ganaderas de las zonas se 
muestran en el Cuadro 1. 



Cuadro 1 Parametros productivos en fincas ganaderas (n=24) de Santo Domingo, Ecuador. 2006. 



Parametro 


Promedio (±Se) 


Carga Animal (UA/ha) 


10.63 ±8.69 


Produccion de leche (l/dia) 


146.28 ± 29.14 


Desparasitacion interna (veces al ano) 


2.91 ± 0.26 


Desparasitacion externa (veces al 
ano) 


17.79 ± 1.98 


Concentrado (kg/animal/dia) 


1 .6 ± 0.1 


Banano (kg/animal/dia) 


2.6 ±0.5 


Melaza (l/animal/dia) 


0.2 ± 0.07 


Sal mineral (kg/animal/dia) 


0.09 ±0.9 


Sal comun (kg/animal/dia) 


0.04 ±0.01 



El numero de animales y la composicion del hato fue muy variable en las fincas ganaderas de 
la zona, con fluctuaciones entre 1 5 y 540 cabezas de ganado por finca (Cuadro 2). 




Revista 

CIENCIA 



Vol. 12, 1, 67-80 (2?q9J 



Cuadro 2 Numero y tipo de animales presentes en tineas ganaderas (n=24) de Santo Domingo, 

Ecuador. 2006. 



Tipo de ganado 


Promedio 

(+Se)/finca 


Rango 


Vacas en produccion 


24.9 ± 1.7 


0-90 


Vacas secas 


23.9 ±0.6 


0-55 


Terneras amamantando 


22.8 ±0.9 


0-90 


Reemplazos 1 (destete-lera monta) 


7.4 ± 1.1 


0-15 


Reemplazos 2 (1 era monta- 1 er parto) 


14.5 ± 0.6 


0-45 


Animales de engorde (menores de 1 ano) 


71 .88 ± 0.8 


0-207 


Toros 


33.6 ± 1.1 


0-540 



En las tineas ganaderas de Santo Domingo, apenas el 16.7 % de los productores han 
realizado cruces entre animales, con la finalidad de encontrar una cruza que se adapte mejor a 
las condiciones climatologicas de la zona y que mantenga su potencial productivo. Es asi que 
se puede encontrar tipos de cruzas como por ejemplo Brahman + Criollo, Brown swiss + 
Criollo, Brown swiss + Mestizo, Brown swiss + Holstein (Cuadro 3). 

Cuadro 3 Frecuencia de razas y cruzas de ganado presentes en tineas ganaderas de Santo Domingo, 

Ecuador. 2006. 



Razas y cruces 


No. Fincas con la raza 


Brahman 


6 


Brahman + Criollo 


1 


Brown swiss 


4 


Brown swiss + Criollo 


1 


Brown swiss + Mestizo 


1 


Brown swiss + Holstein 


1 


Cebu 


1 


Criollo 


5 


Gir lechero 


1 


Girolando 


2 


Holstein 


3 


Jersey 


1 


Mestizo 


3 



En la zona de estudio el 94.4 % de las tineas ganaderas poseen pasto Brachiaria. El 
alto porcentaje de productores que prefieren pasto Brachiaria debe posiblemente a que 
representa un cultivo economico ya que no requiere de labores agricolas como fertilizacion y 
control de plagas. La unica labor que se realiza en este pasto es el control de malezas que 
puede ser manual (practicada por el 29.16 % de los productores a una frecuencia que varia 
entre 1 y 12 veces por ano) o quimico (realizada por el 8.33% de los finqueros, bajo una dosis 
promedio de 0.5± 0.9 l/ha). 

Ademas, los productores tambien tienen pastos mejorados en sus potreros como 
saboya, (presente en el 58.33 % de las tineas encuestadas), gramalote y pasto elefante (en el 
12.50 %), pasto miel (en el 8.33 %) y pasto estrella, king grass y mala alfalfa, presentes en el 
4.17 % de las tineas. El 8.33% de los productores aplican fertilizante quimico en los pastos 
mejorados a una dosis promedio de 45 (kg/ha). 

3.1.3 Componente arboreo 

La tineas ganaderas de Santo Domingo de los Colorados tienen arreglos forestales de 
diferente indole. El arreglo mas abundante constituyen los hueros frutales que se encuentran 
presentes en el 70.83 % de las tineas, seguido de los bosques secundarios presentes en el 
41.47 % de las tineas, los bosques primarios y riparios en el 45.1 % de las tineas y las 
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plantaciones forestales en el 33.33 % de las fincas El tamano promedio de los arreglos 
presentes en las fincas ganaderas se puede observar en el Cuadro 4. 



Cuadro 4 Area, error estandar y rango de los arreglos presentes en fincas ganaderas de Santo 
Domingo, Ecuador. 2006 (De acuerdo a los productores). 



Arreglo 


Area promedio por finca (ha) 


Rango (ha) 


Potrero 


76 + 5.26 


2.5-109 


Pasto mejorado 


30.52 + 7.20 


0.3-5 


Pasto natural 


42.99 + 7.77 


1-109 


Pasto code 


0.96 ±0.47 


0.35-0.5 


Cultivo 


11.71 ±4.29 


1 -50 


Bosque primario 


15.4416.40 


0.01 -49 


Bosque secundario 


6.99 ±3.32 


0.5-5 


Plantaciones forestales 


2.111 1.01 


0.2-75 



Las especies mas comunes que se encontraron en los bosques primarios, secundarios 
y riparios fueron el pachaco ( Schizolubium parahybum), el laurel (Cordia alliodora) y la teca 
(Tectona grandis). En menor proporcion se encontraron especies como el Tangare (Carapa 
guianensis), Copal ( Bursera cuneata ), y Moral (Chlorophora tinctoria). 

3.1 .3.1 Arboles en potreros 

En la zona de estudio se encontro que en la mayorfa de fincas dejan arboles en 
potreros (91.67 %). El 40.9 % de los productores senalaron que los arboles presentes en 
potreros fueron de regeneracion natural y el 59.1 % de los productores mencionaron que han 
sembrado arboles en potreros. La principal especie presente en los potreros fue la guayaba, 
utilizada principalmente como fruto para consumo y tambien para el ganado. Ademas existe un 
alto porcentaje de limon y naranja que son utilizados como frutales para consumo y sombra 
para el ganado. 

Los productores ganaderos de la zona reconocieron el valor de los arboles en sus 
fincas y los beneficios de mantener arboles en potreros ya que constituyen fuente de postes, 
madera y lena para la finca y brindan servicios como sombra y proveer frutos a los animales. La 
mayor cantidad de arboles en potreros reportados en la zona fueron usados como frutales (por 
el 100 % de los productores), los cuales fueron generalmente utilizados para autoconsumo. 
Dentro de estas especies representatives de este grupo se encontraron: Limon (Citrus 
reticulata), presente en el 62.5 % de las fincas, Naranja ( Citrus sinensis, en el 58.3 %), 
Guayaba ( Psidium guajaba, en el 54.2 %), y Guanabana ( Annona muricata, en el 37.5 %). 

En segundo orden se encontraron los arboles que se emplean como sombra para el 
ganado (por el 83.33 % de los productores). Luego se encontraron los arboles que fueron 
utilizados como madera, por el 75 % de los finqueros. En menor escala se ubicaron los arboles 
usados como postes muertos (Figura 1). 
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Figura 1 Frecuencia del uso de los arboles en potreros, Santo Domingo, Ecuador. 2006. 

Las especies de arboles maderables mas comunes en potreros fueron el Mamey y 
Laurel presentes en el 62.5 % y 54.17 % de las tineas encuestadas, respectivamente. Otra 
especie maderable frecuente en la zona fue el Guayacan (presente en el 41 .67 % de las 
tineas). 



Otro uso importante de los arboles en las tineas fue como sombra para el ganado; los 
arboles comunmente utilizados para este proposito fueron los frutales, como el limon (presente 
en el 41.7 % de las tineas), la naranja y la guayaba (presentes en el 25 % de las tineas). 

3.1 .3.2 Cercas vivas 

Se encontro que en el 100 % de las tineas tuvieron cercas vivas. El 50 % tuvo cercas 
vivas compuestas por una sola especie y el 50 % de las tineas tuvo cercas vivas compuestas 
por mas de una especie. Segun estimaciones de los productores, la longitud total promedio de 
las cercas vivas por tinea fue de 1 1 .44 ± 1 5.84 km con variaciones entre 0 y 70 km. 

Las principales especies arboreas registradas en cercas vivas fueron Caraca (Erythrina 
poeppigiana) presente en el 62.5 % de las tineas y Nacedero ( Trichanthera gigantea) en el 
16.67 %. Otras especies que se encontraron en menor frecuencia fueron el Laurel y Yuca de 
raton presente en el 2.8 % de las tineas. 

Los arboles presentes en las cercas vivas generalmente fueron podados con una 
frecuencia de una y 12 veces por ano. El 44.67% de los productores podaron los arboles de las 
cercas una vez al ano, el 25 % lo hicieron dos veces, el 4.17 % tres y el 4.17 % doce veces. 
Los residuos de la poda de caraca fueron generalmente utilizados como alimento para el 
ganado y como material vegetativo para propagacion. 

3.1.4 Componente herbaceo 



En la zona el principal pasto encontrado en los sistemas ganaderos fue la brachiaria 
presente en el 91 .67 % de las tineas y cubriendo el 45.72 % del area total de potreros. Tambien 
se encontro saboya (P. maximun), presente en el 62,5% de las tineas. Cabe anotar que 
saboya es el pasto mas difundido en el litoral, sin embargo en las 24 tineas registradas el pasto 
predominante fue Brachiaria (Cuadro 5). 
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Cuadro 5 Especies de pastos dentro de las fincas, Santo Domingo, Ecuador, 2006. 



Especie 


Nombre cientifico 


No.de Fincas 


Brachiaria 


Brachiaria sp. 


22 


Elefante 


Pennisetum purpureum 


2 


Estrella 


Cynodon nlemfluensis 


2 


Grama 


Paspalum macrophylum 


1 


Gramalote 


Paspalum fasciculatum 


2 


King grass 


Saccharum sinense 


1 


Maralfalfa 


Pennisetum violaceum 


2 


Manf forraiero 


Arachis pintoi 


2 


Micay 


Axonopus micay 


1 


Miel 


Paspalum dilatatum 


4 


Saboya 


Panicum maximun 


15 



3.2 Caracteristicas de los productores 



Los productores ganaderos de la zona mencionaron haber tenido experiencia en 
actividades ganaderas entre 1 y 50 anos, lo cual muestra la heterogeneidad de los finqueros en 
esta actividad. El tiempo promedio de los productores en la zona de Santo Domingo fue de 
16.4 anos con variaciones entre 1 y 50 anos. En promedio los productores ganaderos utilizaron 
como mano de obra familiar 1.15 personas/finca y 4.48 personas/ finca como mano de obra 
contratada (Cuadro 6). 



Cuadro 6 Utilizacion de mano de obra en las fincas, expresado como numero de personas/ano que 



trabajan en la finca. Santo Domingo de los Colorados, Ecuador, 2006. 



Variable 


Promedio 


MO familiar (personas/ano) 


1 .15 ± 0.41 


MO contratada (personas/ano) 


4.48 ± 1.24 


% familias en las que algun miembro trabaja fuera de la finca 


50 



3.3 Relaciones de las caracteristicas de las fincas con la cobertura arborea 



La cobertura arborea total (CA) (area de bosque primario + area de bosque secundario 
+ area de plantaciones forestales + area de bosque ripario) de las fincas ganaderas en la zona 
de Santo Domingo de los Colorados estuvo explicada por 6 modelos de regresion lineal simple; 
en forma directamente proporcional por el % de area de cultivos (P=0.0380), e inversamente 
proporcional por el area de potreros (P=0.0738), la produccion de leche (P=0.0326), la edad 
del destete (P=0.0189), el tiempo de pastoreo (P= 0.0288) y los anos que los productores se 
han dedicado a la ganaderfa (P= 0.0003; Figura 2). 
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Figura 2 Relacion entre el porcentaje de cobertura arborea total y las variables socioeconomicas. 
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4 DISCUSION 

4.1 Uso del suelo 

En la zona de Santo Domingo de los Colorados predominan las tineas doble proposito 
con una superficie promedio de 93,7 ha. Sin embargo, existen tineas de Carne, leche y tineas 
dedicadas a la agricultura. Las tineas de la zona de Santo Domingo de los Colorados tienen un 
area promedio de 90,96 ha, con variaciones entre 5 y 300 ha. Estas variaciones se pueden 
deber a que los propietarios de tineas ganaderas con grandes extensiones de terreno, 
adquieren propiedades de menor tamano para incrementar sus predios, provocando que los 
productores con predios de menor tamano se encuentren en menor cantidad dentro de la zona. 



4.2 Caracteri'sticas de los productores 

El 54.2 % de los productores de la zona de Santo Domingo viven en la tinea y son los 
encargados de tomar las decisiones en el predio. Solamente el 25 % de los productores tienen 
a la tinea como fuente principal de ingresos y el 75 % se dedican a otra actividad simultanea a 
la produccion agropecuaria, tales como: empleos privados, militares, transportistas, 
comerciantes, entre otros. La mano de obra familiar promedio en las tineas de la zona es de 
1.15 jornales/ano y el 54,2 % de las tineas utiliza algun miembro de la familia en las actividades 
del predio. El 70,8 % de las fincas de la zona utiliza mano de obra contratada, en promedio de 
4.48 jornales/ano. 



4.3 Arboles en potreros 

La presencia de arboles en potreros en Santo Domingo de los Colorados, fue 
encontrado en el 91.67% de las fincas encuestadas. Similares resultados se reportaron en 
otros estudios llevados a cabo en Rfo Frio, Costa Rica, en donde el 97.5 % de las fincas 
presentaron arboles en potreros; de igual manera en otro estudio realizado en La Fortuna, 
Costa Rica se encontro que mas del 90% de las fincas ganaderas presentaron arboles 
dispersos en los potreros (Souza de Abreu et al, 2000). En otro estudio llevado a cabo en la 
region Caribe de Colombia se encontro que entre el 26 y 69% de las pasturas de las tineas 
tienen arboles (Cajas-Giron y Sinclair, 2001). 

La mayor cantidad de arboles en potreros reportados en la zona de Santo Domingo de 
los Colorados son frutales. Estos resultados son diferentes a los encontrados en otro estudio 
desarrollado en fincas ganaderas de Costa Rica en donde se encontro que la mayor cantidad 
de arboles presentes en potreros son maderables (Villafuerte, 1998; Morales y Klienn, 2000). 
En cuanto a la presencia de arboles maderables en potreros se encontro que las especies mas 
abundantes fueron el Mamey y el Laurel. 

La cantidad de especies de arboles en potreros es muy variable entre tineas. En la 
zona de Santo Domingo de los Colorados los productores reportaron 41 especies diferentes de 
arboles diferentes dentro de los potreros. Estos resultados son diferentes a los encontrados en 
otros estudios llevados en otras zonas como Canas, donde se reportaron que 60 especies 
diferentes de arboles estan presentes en potreros (Morales y Kleinn, 2001). En la misma zona, 
en otra investigacion realizada, se reportaron que 64 especies forman parte del componente 
arboreo presente en potreros (Restrepo 2002). 

Los principales usos de los arboles presentes en potreros por los finqueros fueron 
como frutos (95.8 % de los finqueros), como forraje para el ganado (41 .7 %), sombra (29.2 %) 
y en menor escala como lena, madera, postes y forraje para el ganado; estos usos fueron 
anteriormente reportados en fincas ganaderas de otra zona, sin embargo el uso mas frecuente 
reportado fue como madera (Villafuerte, 1998). 
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4.4 Cercas vivas 

En la zona de Santo Domingo de los Colorados, el 100 % de los productores reportaron 
que tienen cercas vivas en sus tineas, con predominancia de la especie Caraca. Similares 
resultados se encontraron en las zonas de Rio Frio y La Fortuna, Costa Rica en donde el 
87.3% y el 85 % de las tineas ganaderas, respectivamente, tuvieron cercas vivas compuestas 
principalmente por la misma especie (Villacis, 2003; Souza de Abreu etal., 2000). 



El 44.67% de los productores podan los arboles de las cercas una vez al ano, 
principalmente los arboles de Caraca y Nacedero. La poda se realiza en cualquier epoca del 
ano independientemente de la presencia de lluvias. Los residuos de la poda son utilizados para 
brindar alimento al ganado y para cortar nuevas estacas usadas como material vegetativo para 
propagacion Ningun productor de la zona reporto que usa las especies arboreas de las cercas 
vivas como fuente de madera o lena, contrario a la zona de Canas donde los productores 
utilizan los arboles de las cercas vivas con fines maderables (Harvey et al., 2003; Restrepo 
2002 ). 



Otras especies que reportaron los finqueros que forman parte de las cercas vivas 
presentes en sus tineas fueron: Laurel, Yuca de raton. Sin embargo estas especies no son 
manejadas o su manejo es minimo en comparacion a la caraca y nacedero. 

4.5 Relaciones de las caracteri'sticas de las fincas con la cobertura arborea 

Al analizar detalladamente cada una de las ecuaciones, el area de potreros es 
inversamente proporcional e influye negativamente en sobre la cobertura arborea, debido a que 
al incrementar las areas de pastoreo en las fincas ganaderas se eliminan los remanentes de 
bosque presentes. El area de cultivos presento una relacion directamente proporcional e influye 
positivamente en la cobertura arborea, esto se podria explicar ya que algunos cultivos como el 
cacao y el cafe requieren de la provision de sombra de los arboles para poder cumplir su ciclo 
productivo. La produccion de leche influyo negativamente en la cobertura arborea, esta relacion 
se podria deber a que al incrementar las areas de pasto mejorado que tienen altos niveles de 
nutrientes se debe eliminar los arboles de los potreros, ya que la sombra ejercen un efecto 
negativo sobre los pastos (Villacis 2003). La edad del destete indico una relacion inversamente 
proporcional e influye negativamente en la cobertura arborea, debido a que conforme se van 
destetando animales, el ganadero puede verse en la necesidad de aumentar el area dedicada 
a la ganaderia, consecuentemente diminuyen las areas boscosas. El tiempo del pastoreo 
sugirio una relacion negativa en la cobertura arborea, ya que el ganado podria estar pisoteando 
los arboles de la regeneracion natural de los potreros. Finalmente, los anos que los propietarios 
se han dedicado a la ganaderia tambien presentaron una relacion negativa con la cobertura 
arborea, esto se puede explicar ya que los sistemas productivos tradicionales han sugerido la 
eliminacion del componente arboreo para la implementacion de potreros. Estos resultados son 
diferentes a los reportados en otra investigacion desarrollada en la zona de Rio Frio, Costa 
Rica en donde la cobertura arborea total de las fincas estuvo explicada por un modelo de 
regresion multiple en funcion del area total de la finca, la mano de obra contratada, el area de 
potreros, la frecuencia de desparasitacion al ganado, el tiempo de ocupacion de los potreros y 
el suministro de concentrado al ganado (Villacis, 2003; Villacis et al., 2003); y a los encontrados 
en otro estudio llevado a cabo en la zona de Canas, donde la cobertura arborea en potreros 
estuvo explicada por el control de malezas, la mano de obra familiar, y el periodo de descanso 
en los potreros (Restrepo, 2002). 

5 CONCLUSIONES 

Los tipos de cobertura arborea predominantes en fincas ganaderas de Santo Domingo de los 
Colorados constituyen los arboles dispersos en potreros presentes en el 95.7 % de las fincas y 
las cercas vivas presentes en el 1 00 % de las fincas ganaderas de la zona. 
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Los productores ganaderos reconocen el valor de los arboles en sus fincas y las ventajas de 
mantener arboles en potreros. Entre las principales ventajas son los productos (que proveen 
como postes, madera, lena, forraje, frutos) y servicios como brindar sombra y proveer frutos a 
los animales. 



Las especies arboreas mas representativas que se encuentran presentes en potreros, 
mencionadas por los productores fueron el limon (Citrus reticulata, presente en el 62.5 % de las 
fincas), la naranja ( Citrus sinensis, en el 58.3 %), la guayaba (Psidium guajaba, en el 54.2 %), y 
la guanabana (Annona muricata, en el 37.5 %). 

Bajo las condiciones de estudio el porcentaje de area de cultivos influyo positivamente sobre la 
cobertura arborea, el area de potreros, la produccion de leche, la edad del destete, el tiempo 
de pastoreo y los anos de los productores dedicados a la ganaderia influyeron en 
negativamente la cobertura arborea total presente en las fincas. 



La zona de Santo Domingo de los Colorados se dedica a la actividad ganadera y representa un 
potencial en la produccion ganadera del Ecuador. Por tradicion los productores han utilizado 
diferentes especies de arboles de sus fincas en las actividades productivas, generalmente para 
la provision de sombra y frutos al ganado y en algunas ocasiones para postes. Sin embargo, 
los arboles no reciben ningun tipo de manejo, con excepcion del componente arboreo presente 
en cercas vivas el cual es manejado por los productores para disminuir la sombra en los pastos 
y para obtener forraje para el ganado. 

En la mayoria de las fincas los trabajos que se puedan proponer en cuanto a manejo e 
implementacion de sistemas silvopastoriles podrian ser satisfactorios siempre y cuando exista 
el apoyo de instituciones financieras que faciliten la provision de arboles y el personal tecnico 
necesario para llevar as cabo estas actividades. Por estas razones, el conocimiento previo de 
una zona es de suma importancia para el desarrollo de proyectos de investigacion que esten 
enfocados a las condiciones y caracteristicas de los productores. La informacion recopilada en 
este trabajo muestra las principales caracteristicas sobre las fincas ganaderas de la zona en 
cuanto a uso del suelo, caracteristicas de los productores, niveles de intensificacion, tipologias 
de los sistemas productivos, presencia y utilizacion del componente arboreo en las fincas. La 
informacion generada en este estudio se deberfa utilizar como base para la implementacion de 
sistemas silvopastoriles de tal forma de rescatar y aumentar el componente arboreo y de esta 
forma evitar que se pierdan los pocos arboles que se encuentran presentes en las fincas 
ganaderas de la zona. 
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SOBRE EL EQUILIBRIO DE ESTRUCTURAS APORTICADAS 



Pablo Caiza Sanchez 
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Escuela Superior Politecnica del Ejercito 

RESUMEN 

Se usa la teoria completa no lineal de vigas planas para analizar un portico piano de 
un piso y un vano, con conexiones rigidas y fijo en su base. Las ecuaciones de equilibrio no 
lineales se plantean mediante principios variacionales usando el metodo de elementos finitos 
con polinomios lagrangianos de diferente orden como aproximacion de los campos 
cinematicos. El metodo de interpolacion de Newton con control de longitud de arco fue usado 
para la solucion de estas ecuaciones no lineales. Diferentes tipos de cargas asf como 
dimensiones geometricas fueron implementados con el fin de evaluar los resultados de 
desplazamientos y fuerzas internas. El efecto de cambiar el modulo elastico del material 
tambien fue considerado. Finalmente se tomaron en cuenta los efectos de imperfecciones 
geometricas. 



ABSTRACT 

A plane one-bay-one-story frame with rigid connections and fixed at its base is 
analyzed using the fully nonlinear planar beam theory. The nonlinear equilibrium equations 
were written using variational principles and the finite element method. The kinematic fields 
were approximated by Lagrangian polynomials. The Newton method with arc-length control 
was used for the solution of the nonlinear equilibrium equations. Different loading scenarios as 
well as geometric dimensions were implemented in order to evaluate displacement and internal 
force results. The effect of changing the material moduli was also considered. Finally geometric 
imperfections were also taken into account. 

1 INTRODUCCION 

Este articulo trata de la metodologia para resolver estructuras aporticadas usando la 
teoria no lineal completa de vigas planas. Es importante recalcar que la no linearidad se 
refiere a la ecuaciones de equilibrio y deformacion-desplazamiento y no a las 
ecuaciones constitutivas. Cinco pasos basicos son necesarios: 

Descripcion de la geometria del sistema estructural. 

Descripcion de las cargas externas aplicadas. 

Establecimiento de las ecuaciones de equilibrio en la configuracion deformada. 

Aplicacion del metodo de elementos finitos. 

Uso del metodo de interpolacion de Newton con control de longitud de arco. 

Esta metodologia se aplico a la solucion de un portico piano de un piso y un vano con 
diferentes caracteristicas geometricas asf como de carga. Para probar la validez de los 
resultados algunas ecuaciones de equilibrio son descritas con mas detalle. Adicionalmente los 
resultados son comparados con aquellos obtenidos con el programa comercial SAP2000. 

2 DESCRIPCION DEL SISTEMA ESTRUCTURAL 

La estructura basica a ser analizada es descrita en la Figura 1 que se indica a 
continuacion. 
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Figura 1 Estructura basica a ser analizada. Fuente: Handout CEE471 Term Project 2 [2] 

Dos casos son implementados, el primero llamado “Columnas Esbeltas” y el otro “Viga 
de Gran Longitud”. Las dimensiones y propiedades geometricas mas importantes, asi como el 
valor de los modulos de elasticidad y de corte usados, se encuentran en la Tabla 1 . 



Tabla 1 Caracterfsticas de los dos casos implementados. 



E=2 10 000T/m2 


G=1 50 000 T/m2 




“Columnas Esbeltas” 


“Viga de Gran Longitud” 


H (m) 


7.0 


3.5 


L(m) 


4.5 


9.0 


Ic (m4) 


0.000675 


Ac (m2) 


0.09 


lb (m4) 


0.000893 


Ab (m2) 


0.0875 



Se observa que los valores de E y G son una decima parte de los valores para 
hormigon armado. Lo anterior se usa para producir mayores deformaciones y hacer mas 
evidente el comportamiento no lineal de la estructura. Notese ademas que H es la altura de las 
columnas y L la longitud de la viga. Las cargas aplicadas son de dos tipos: vertical y lateral. 
Cada uno de estos tipos corresponde a un caso de carga, siendo un tercer caso, la 
combinacion de ambos tipos. Los tipos de carga, en el paso inicial (cero) de analisis, son 
descritos en la Tabla 2. 



Tabla 2 Casos de carga 



Caso Vertical 


qbl 


9 T/m (constante) 


qb2 


0 T/m (variable) 


Pv col 


2 T (variable) 


Caso Horizontal 


Ph col 


1 T (variable) 



Se observa que hay cargas constantes qbl y otras variables qb2 (carga distribuida 
sobre la viga), Pvcol (carga vertical puntual sobre las columnas), Phcol (carga horizontal 
puntual sobre las columnas). Las cargas variables se emplean para determinar nuevas 
configuraciones de equilibrio bajo cargas diferentes. 

3 PROCEDIMIENTO ANALITICO 

3.1 La teoria no lineal completa para vigas planas 

Las ecuaciones diferenciales que gobiernan el equilibrio en la estructura deformada 

son [1]: 
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H' + p = 0 V' + o = 0 

r ‘ ( 1 ) 
W + V{1 + u'> - Hw' + m=C 

Donde H,V son las fuerzas resultantes axial y transversal al eje de cada elemento 
estructural no deformado, M es el momento resultante; p,q, fuerzas distribuidas transversal y 
axial y m momento distribuido; u es el desplazamiento del centroide en la direccion axial y w el 
desplazamiento del centroide en la direccion transversal. 

Estas ecuaciones son de equilibrio en la estructura deformada. La deformada 
geometrica y el equilibrio de un segmento de viga se muestran en la Figura 2. 




Figura 2 Deformada y equilibrio de un segmento de viga deformada (1 ) 



Usando esta nomenclatura la derivacion de la ecuacion de equilibrio de momentos, por 
ejemplo, es la siguiente: 

M (x + Ax) - M (a') +V(x + A\')(Av + u (x + Ax) - u (x)) + 

- H (x + Ax)(w(x + Ax) - w(x)) + 

Ax 

+ J[(w(£)- w{x))p{£)~ (£ + «(£))?(£) + m{g)]d£; = 0 
0 

Donde ^ mide distancia (en la configuracion no deformada) desde el extremo izquierdo 
del segmento. Dividiendo la ecuacion para Ax da: 
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M (x + Ax) - M (.x) V(x + Ax)(Ax + u(x + Ax) - u(x)) 

Ax Ax 

H(x + Ax)( w{x + Ax) - vv(.x)) 

Ax 

Ax 

+ Y± x J[M#)- w{x))p{^)-{^ + «(£))<?(£)+ = 0 

0 

Tomando el limite Ax— >0, y reconociendo la definicion de derivada se tiene 
nuevamente la ecuacion de equilibrio de momentos que es la tercera de las ecuaciones (1). 

Las fuerzas en una seccion transversal pueden tambien ser expresadas como componentes 
relativas a la seccion transversal rotada: fuerza axial N(x) y fuerza cortante Q(x) [1], 



H = Ncosd— QsinQ 
V = NsinS + Qcosd 



Siendo 0 es la rotacion de la normal al piano de la seccion transversal. 

Estas ecuaciones diferenciales de equilibrio se resolveran mediante principios 
variacionales y el metodo de elementos finitos. Notese por tanto que para cada elemento 
estructural la ecuacion que se desea solucionar es la ecuacion variacional. 

G{A,u,u) = 0 , \/ueT^(0,l) ( 3 ) 

Donde A es un parametro de carga y permite usar cargas tales como q(A,x)=q1 (x)+ A 
q2(x); u = {u(x),w(x),0(x)}] r,(0,/) es un conjunto de funciones m={w,w\#} que 

satisfacen condiciones de borde esenciales y cuyas primeras derivadas son integrables de 
segundo grado. 

De acuerdo a lo anterior se tiene: 

G(A,u,u) = J (EIZ 0 Z 0 + GA P () P {) + EA£ 0 £q \x-W e {X) ( 4 ) 

Donde W / I7 (/l)es el funcional por trabajo externo de las cargas aplicadas. Las 

deformaciones unitarias reales , {3 0 , £ 0 son funciones de los desplazamientos reales, los 
cuales, a su vez, son funciones de la variable independiente x. Similarmente las deformaciones 
unitarias virtuales % 0 ,J3 0 ,£ 0 son funciones de los desplazamientos reales y virtuales, los 
cuales, a su vez, son funciones de la variable x. 

Las deformaciones unitarias son: 



fie = w r cost ? — (1 + u r )sm6 
E a = w'sirtd + (1 + i £,'}cos8 — 1 



Donde xo (la curvatura) ,(3o y Go son las deformaciones unitarias asociadas con M, Q 
y N, respectivamente. 
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Se considera que la viga es elastica y que las relaciones constitutivas son lineales por 

lo tanto: 



M = El'xo Q = GAfio N = EAeo 



( 6 ) 



Siendo El, GA y EA son las constantes de proporcionalidad tradicionales. 



3.2 Representation de los campos de desplazamiento 

Los campos de desplazamiento u, w y 0 son representados por medio de funciones de 
forma lagrangianas. 

Una funcion lagrangiana general de orden n que pasa a traves de cero en n-1 puntos 
especificados x, y tiene un valor unitario en un punto especificado x k tiene la siguiente ecuacion 
[1]: 

x -x. 

= U*k x -x. (7) 

k i 



l ( x) = n tl 

n i 



Como un ejemplo las funciones de forma cuadraticas (correspondientes a polinomios 
lagrangianos de segundo grado) son descritas en la Figura 3 a continuacion: 




Figura 3 Funciones lagrangianas de interpolacion de segundo grado. 
Fuente: Handout CEE471 Term Project 2 [2] 



Observese que aqui u(x) es una funcion general que puede representar 
desplazamientos axiales, transversales e inclusive giros. Ademas notese que el elemento 
estructural (viga o columna) necesita en este caso tres nudos para acomodar estas funciones 
cuadraticas. 

3.3 El metodo de interpolacion de Newton con control de longitud de arco 

Cada configuracion de equilibrio se determino usando el metodo de interpolacion de 
Newton con control de longitud de arco. Este se empleo repetidamente cuando las condiciones 
de carga se hicieron variar. 

Algunas de la mas importantes caracteristicas del metodo de interpolacion de Newton 
aplicado a este problema son descritas en los siguientes parrafos: 

Definense las matrices A(A,a) y b(A,a) de la siguiente manera [1]: 



Ai;A, a) = 



K{a) 

.2 (a - a n ) 



fefT> 

2(A — /L w 
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£>( 2 , a) = 



-g(A, o) 

-r(l, a) 



Donde a es el vector que contiene los coeficientes de las funciones de interpolacion 
que se calculan en la iteracion v del caso de carga n+1 del proceso de Newton, a n es el vector 
que contiene los coeficientes de las funciones de interpolacion calculados en el caso de carga 
previo, An el factor de carga del paso previo, A el factor de carga empleado en el paso en 
consideracion, K(a) es la matriz de rigidez tangente, k(2) es el vector conteniendo las cargas 
distribuidas y puntuales. Notese que b(Jl,a) es el vector que contiene la cantidad por la cual el 
equilibrio y la restriccion de la longitud de arco fallan en satisfacer al estado (A, a). “De alii que, 
cualquier configuracion que no esta a una distancia fija del punto previo y no esta en 
equilibrio puede ser mejorada resolviendo 



A(A,a v ) Ax v =b(A,a v ) 

Siendo Ax=(a,A) es una matriz que guarda la variacion de los coeficientes del estado 
de desplazamiento y el factor de carga empleado. 

La configuracion mejorada puede ser obtenida como [1]: 

x v+1 = x v +Ax v 

“El numero de iteraciones es el necesario para que la solucion este dentro de una 
tolerancia especificada”[1]. 

4 VERIFICACION Y PRECISION COMPUTACIONAL 

La teorfa expuesta previamente fue implementada en un programa en MATLAB. Como 
paso inicial y dado que el metodo de elementos finitos provee una aproximacion de la solucion 
exacta dada por las ecuaciones diferenciales de equilibrio es esencial una investigacion de la 
precision de los algoritmos usados. Para el efecto se emplearon diferentes grados para las 
aproximaciones polinomiales lagrangianas. Notese sin embargo que, puesto que la solucion 
exacta para las funciones de desplazamiento son polinomios de cuarto grado, funciones de 
forma de cuarto orden son suficientes para capturar una solucion precisa. 

Adicionalmente los resultados obtenidos se compararon con los correspondientes a los 
calculados con SAP2000.En la comparacion que se presenta se usaron ecuaciones 
lagrangianas de quinto orden. 

Los resultados mostrados corresponden a la estructura denominada “Columnas 
Esbeltas” unicamente con una carga lateral de 5 T. La Figura 4 muestra la configuracion 
deformada. 
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Figura 4 Configuration deformada usada para validacion. 



Los desplazamientos globales en los nudos superiores del portico fueron chequeados 
contra los resultados obtenidos con el programa comercial SAP2000 (opcion: “large 
displacements”) y se muestran en la Tabla 3. 

Tabla 3 Desplazamientos globales en los nudos en la parte superior del portico ana lizado 





Nodo 2 


Nodo 3 




Viga no lineal 


SAP2000 


Viga no lineal 


SAP2000 


Ax [m] 


0.596 


0.621 


0.596 


0.621 


Ay [m] 


-0.028 


-0.026 


-0.030 


-0.029 


0 [rad] 


0.00335 


0.00332 


0.00323 


0.00332 



Los resultados son similares dando confianza para estimar otros modelos 
estructurales. 

Es importante indicar adicionalmente que funciones de forma de mayor grado dieron 
resultados igualmente similares pero que, en el caso de funciones de forma menores al cuarto 
grado, la estructura se “rigidizo”. Este es un fenomeno numerico conocido como “shear 
locking” debido a que las funciones de forma no son capaces de representar las 
deformaciones exigidas por las cargas sobre la geometria estructural. 

Finalmente otro efecto no lineal importante es el pandeo. De manera similar a los 
resultados de desplazamientos, las cargas de pandeo fueron mas parecidas a la teoricas 
mientras mayor era el grado de las funciones de forma empleadas, siendo sin embargo las 
funciones de cuarto orden suficientemente precisas. La carga de pandeo se determina cuando 
uno de los valores propios de la matriz de rigidez tangente alcanza el valor de cero. Notese 
que este ultimo hecho es posible debido a la teorfa no lineal usada: basicamente las 
ecuaciones de equilibrio se plantean sobre la estructura deformada. 



5 GEOMETRIA Y CARGAS 

5.1 Caso Nol : Columnas esbeltas unicamente con carga vertical 



Este ejemplo corresponde a la geometria denominada “Columnas Esbeltas” y al caso 
de carga “Vertical”. Diez pasos de carga son usados con un factor de magnification de 10, lo 
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que significa una carga final distribuida sobre la viga qb=9T/m y una carga puntual sobre cada 
columna Pvcol=200 T (ver Figura 5). 




Ademas y como un ejemplo de las curvas desplazamiento transversal-cargas, la 
Figura 6 muestra la que se obtuvo para el nudo en el centro de la altura de la columna a la 
izquierda de la estructura. 

Se observa un comportamiento inicial lineal en un rango de desplazamientos pequeno 
y luego, un brusco aumento para de nuevo estabilizarse con un comportamiento 
aproximadamente lineal. Este comportamiento se puede explicar matematicamente 
observando una continua reduccion de los valores de los vectores propios de la matriz de 
rigidez tangente a medida que la carga aumenta hasta llegar a cero, es decir pandeo, y un 
brusco aumento de la carga en este momento. Signo de una estabilizacion debida a que los 
otros elementos de la estructura vienen en la ayuda de la columna. Es un comportamiento 
caracteristico de estructuras hiperestaticas. 




Figura 6 Curva desplazamiento transversal-carga puntual vertical (nudo en H/2, columna a la izquierda) 
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5.2Caso 2: Columnas Esbeltas con cargas verticales y laterales 

Este ejemplo corresponde a la geometna denominada “Columnas Esbeltas" y al caso 
de carga “Vertical” anadido al “Horizontal”. Ocho pasos de carga fueron usados con un factor 
de incremento de 2, lo que significa una carga final qb=9T/m, Pvcol=32 T y Phcol=16 T (ver 
Figura 7). 




Figura 7 Configuracion deformada para el Caso No 2 



El algoritmo implementado se vuelve inestable mas alia de esta maxima carga. 

5.3Caso 3: Viga de gran longitud con cargas verticales y horizontales 

Este ejemplo corresponde a la geometna denominada “Viga de Gran Longitud” y al 
caso de carga “Vertical” anadido al “Horizontal”. Diez pasos de carga son usados con un factor 
de magnificacion de 2, que significa una carga final qb=9T/m, Pvcol=40 T y Phcol=20 T (ver 
Figura 8). 
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Figura 8 Configuracion deformada para el Caso No3 

5.4Caso 4: Cambio del modulo de elasticidad 

Los Casos 2 y 3 fueron nuevamente solucionados pero en este caso con la mitad del 
valor de su modulo de elasticidad (ver Figura 9). 

En ambos casos el algoritmo implementado con los mismos parametros se volvio 
inestable. Por lo tanto por ejemplo en el Caso 2 unicamente diez pasos de carga fueron 
usados con un factor de magnificacion de 1, conduciendo a una carga final de qb=.1T/m, 
Pvcol=2 T y Phcol=1 T. Similarmente se debieron reducir las cargas del Caso 3. 

6 IMPERFECCIONES GEOMETRICAS 

La viga y la columna a la derecha fueron calculados como elementos geometricamente 
perfectos, pero la columna a la izquierda tiene desplazamiento laterales lineales con un maximo 
de 0.1 5 m a H/2 y 0 en sus extremos (ver Figura 1 0). 

A pesar de que la deformada general es similar a anteriores, hay una gran reduccion 
de las cargas para un determinado desplazamiento o un gran aumento de los desplazamientos 
para una carga determinada. 
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Figura 9 Configuracion deformada para el Caso2 con la mitad de su modulo elastico inicial 




Figura 10 Imperfecciones geometricas: la columna a la izquierda esta desplazada 



7 DISCUSION 

Las columnas en el Caso 1 (columnas esbeltas bajo carga vertical) presentaron 
deformaciones extraordinarias las cuales denotan pandeo local seguido por la estabilizacion de 
la estructura. La viga, en este caso, fue mucho menos importante y se puede considerar que 
simplemente “acompano” los desplazamientos de las columnas. 
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Las cargas laterales anadieron nuevas caracteristicas al comportamiento estructural 
del Caso 2 comparado con el caso previo. Los desplazamientos laterales fueron mucho 
mayores y, lo mas importante, las cargas verticales pudieron alcanzar unicamente una fraccion 
del Caso 1 . La diferencia se produce por el efecto P-delta en las columnas. 

Las columnas en el Caso 3 fueron mas cortas que en el Caso 1 o 2. Por otro lado la 
longitud de la viga es mucho mayor. El comportamiento estructural en este caso es mas 
estable, mostrando la importancia de tener columnas menos esbeltas. Sin embargo las 
deflecciones de las vigas fueron 3 veces mayores que las de las columnas. 

Los Casos discutidos anteriormente fueron caracterizados por cargas cerca de la 
razon de convergencia del algoritmo de Newton implementado sobre todo con los parametros 
que definen el control por longitud de arco. De aquf que debilitar al material produjo 
inestabilidad computacional. Sin embargo queda claro que las cargas que soportarfa la 
estructura son mucho mas pequenas y que los problemas de pandeo son mas severos. 

Finalmente, debe notarse que los campos de desplazamiento en el portico con 
irregularidades geometricas fueron mucho mayores que sin imperfecciones geometricas. 

8 CONCLUSIONES 

Cinco pasos basicos fueron necesarios para implementar la teorfa nolineal para vigas 
planas: descripcion de la geometrfa del sistema estructural, descripcion de las cargas externas 
aplicadas, las ecuaciones de equilibrio en la configuracion deformada, aplicacion del metodo de 
los elementos finites y el uso del metodo de Newton con control de longitud de arco. 

Estos pasos fueron aplicados para solucionar una estructura aporticada de un piso y un 
vano con diferentes condiciones geometricas y de carga. 

Las caracteristicas consideradas mostraron la estabilidad de las columnas en un 
portico y su comportamiento de pandeo bajo cargas verticales y horizontales. Este 
comportamiento fue caracterizado por desplazamientos perpendiculares a las cargas axiales 
aplicadas. En el caso en que algunas de las cargas produjeron desplazamientos laterales 
adicionales (cargas horizontales a nivel de la viga), las cargas externas resistidas fueron mucho 
menores que de otra manera porque las fuerzas verticales multiplicadas por los 
desplazamientos laterales produjeron nuevos momentos de segundo orden. El comportamiento 
de la viga fue tambien importante porque, si no por su resistencia, sus deflecciones modificaron 
la respuesta del portico. 

Es importante resaltar que estos resultados son validados por el conocimiento de que 
el software usado fue adecuado. Un buen programa computacional debe tener necesariamente 
alguna alternativa de calculo para comparar la precision de los resultados y su uso debe 
basarse en una buena comprension teorica. 

La aplicacion del conocimiento de la teorfa no lineal debe llevar a conclusiones basicas 
tales como que una fabricacion y construccion precisas deben ser enfatizadas para evitar las 
imperfecciones geometricas en las estructuras, pues estas afectan de manera sustancial el 
comportamiento esructural; pero aun mas importante, permitira el desarrollo de soluciones 
novedosas, que pueden darse con el cada vez mayor uso de nuevos materiales de 
construccion, asf como comprender el comportamiento estructural ante situaciones extremas 
tales como sismos, incendios o ataques terroristas. 

% 1 2 3 4 5 6 7 

% * * 

% | Program NonlinearBeam | 

% j Fundamentals of Structural Mechanics, 2nd Edition | 

% j K. D. Hjelmstad, July 1, 2004 | 

% j Modified by P. Caiza, Dec 7, 2007 | 

% * * 
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%.. Set problem parameters, 

clear; tol = l.e-8; alpha = 10; maxsteps = 10; maxit = 100; 
nnstep = 10; npts = 21; 

fname=input('Please type the name of input file [input.dat]: 

%... Default filename will be input.dat' 
if isempty(fname) 
fname='input.dat'; 
end 

%... fid will have the file id. 'r'for read-only 

fid=fopen(fname,'r'); 

nelem=str2num(fgetl(fid)); 

nnodes=str2num(fgetl(fid)); 

nbasis=str2num(fgetl(fid)); 

ndof=str2num(fgetl(fid)); 

ndm = ndof + 1; XX=zeros(nelem,1); YY=zeros(nelem,1);T=zeros(3*nbasis,3*nbasis); 
indice=zeros(nelem,3*nbasis); loadsT=zeros(3*nelem,3);DT=zeros(nelem,3); 
A=zeros(ndm,ndm); b=zeros(ndm,1); 

for nn=1 :nelem %i-node of each element 

temp=str2num(fgetl(fid)); 

XX(nn)=temp(1); 

YY(nn)=temp(2); 

end 



for nel=1 :nelem 

indice(nel,:)=str2num(fgetl(fid)); 

loadsT(3*nel-2,:)=str2num(fgetl(fid)); 

loadsT(3*nel-1,:)=str2num(fgetl(fid)); 

loadsT(3*nel,:)=str2num(fgetl(fid)); 

phiT(nel)=str2num(fgetl(fid)); 

DT(nel,:)=str2num(fgetl(fid)); 

xlengthT(nel)=str2num(fgetl(fid)); 

end 



% numeration for unknowns 
%p1, p2, H(L) 

% ql, q2, V(L) 

% ml, m2, M(L) 

% Transformation angle phi 
% Element moduli D = [El, GA, EA] 
% Element length is xlength 



%.. Initialize values for load step zero, set next trial state 
x = zeros(ndm,1); xo = zeros(ndm,1); x(ndm) = alpha; 
xp=zeros(3*nbasis+1,1); xop=zeros(3*nbasis+1,1); 

%.. Compute MAXSTEPS points along the Equilibrium Path 
for n = 1 :maxsteps 

%.... Perform Newton iteration at each load step 
nu = 0; test = 1.0; 
while (test > tol) & (nu < maxit) 
nu = nu + 1 ; 

% Compute Hessian and residual 

for nel=1 :nelem 
D=DT(nel,:); 

loads=loadsT(3*nel-2:3*nel,:); 

xlength=xlengthT(nel); 

phi=phiT(nel); 

%.. Transformation Matrix 
To=[cos(phi) sin(phi) 0; 

-sin(phi) cos(phi) 0; 
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0 0 1]; 
for i=1 :nbasis 

T(3*i-2:3*i,3*i-2:3*i)=To; %complete T 
end 

for i=1 :3*nbasis 
if indice(nel,i) ~= 0 
xp(i)=x(indice(nel,i)); 
xop(i)=xo(indice(nel,i)); 
else 
xp(i)=0; 
xop(i)=0; 
end 
end 

xp(3*nbasis+1)=x(ndm); 

xop(3*nbasis+1)=xo(ndm); 

[Ap,bp] = fcn(D,xp,xop, loads, xlength,nbasis, alpha, T); 
for i=1 :nbasis*3 
for j=1:nbasis*3 

if (indice(nel,i)~= 0) & (indice(nel,j)~=0) 

A(indice(nel,i),indice(nel,j))=A(indice(nel,i),indice(nel,j))+Ap(i,j); 

end 

end 

if indice(nel,i)~= 0 

b(indice(nel,i))=b(indice(nel,i))+bp(i); 

A(indice(nel,i),ndm)=A(indice(nel,i),ndm)+Ap(i,3*nbasis+1); 

A(ndm,indice(nel,i))=A(ndm,indice(nel,i))+Ap(3*nbasis+1,i); 

end 

end 

A(ndm,ndm)=A(ndm,ndm)+Ap(3*nbasis+1,3*nbasis+1); 

b(ndm)=b(ndm)+bp(3*nbasis+1); 

end 

% Compute residual norm, min eigenvalue of tangent stiffness 

test = norm(b); 

% Compute increment and update state vector 

dx = -A\b; x = x + dx; 

end % while for Newton loop 

%.... Output results, set converged state, guess at next state 

results(x, indice, nelem, XX, YY,xlengthT,n,nu, test, ndm,nbasis,nnstep,npts,phiT) 
temp = xo; xo = x; x = 2*x - temp; 

end % loop on n Load Steps 

%.. End of program NonlinearBeam 



% 1 2 3 4 5 6 7 

% Compute A and b for NonlinearBeam element FCN 

function [A,b] = fcn(D,x,xo, loads, xlength,nbasis, alpha, T) 
ndm=3*nbasis+1; 

%.. Intialize A and b to zero 

b = zeros(ndm,1); A = zeros(ndm,ndm); z = 0; 
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%.. Set weights for Simpson integration 
wt = [1,4,2,4,2,4,2,4,2,4,2,4,2,4,2,4,2,4,2,4,1]; 
intpts = length(wt); dz = 1 /(intpts-1 ); 

%.. Transform x to local system 
x(1 :ndm-1)=T*x(1 :ndm-1); 
xo(1 :ndm-1)=T*xo(1 :ndm-1); 

%.. Loop on integration points 
for n = 1 :intpts 
factor = wt(n)*dz*xlength/3; 

%.... Compute displacements and derivatives 
du = 0; dw = 0; dtheta = 0; theta = 0; 
for i = 1 :nbasis 

[h,dh] = basis(i,z,xlength,nbasis); 
du = du + x(3*i-2)*dh; 
dw = dw + x(3*i-1)*dh; 
dtheta = dtheta + x(3*i)*dh; 
theta = theta +x(3*i)*h; 
end 

ct = cos(theta); st = sin(theta); 

%.... Compute axial strain, shear strain, and curvature 
epsi = dw*st + (1 +du)*ct - 1 ; 
beta = dw*ct - (1+du)*st; 
curv = dtheta; 

%.... Compute axial force, shear force, bending moment, and other forces 
bend = D(1)*curv; shear = D(2)*beta; axial = D(3)*epsi; 

Hor = axial*ct - shear*st; Ver = axial*st + shear*ct; 

Xi = (1+du)*Hor + dw*Ver; Yi = dw*Hor - (1+du)*Ver; 

force = [Hor; Ver; bend]; fyi = [0; 0; Yi]; 

cl = loads(:,1) + loads(:,2)*x(ndm); 

c2 = fyi - cl ; 

c3 = loads(:,2); 

c4 = loads(:,3); 

%.... Compute components of [E(T)DE+G] store in matrix G 



G(1 ,1 ) 


= D(2)*st*st + D(3)*ct*ct; 




G(1 ,2) 


= ct*st*(D(3)-D(2)); 






G(1,3) 


= D(2)*st*(1+epsi) + 


D(3)*ct*beta 


-Ver; 


G(1,4) 


= 0; 






G(2,1) 


= G(1 ,2); 






G(2,2) 


= D(2)*ct*ct + D(3)*st*st; 




G(2,3) 


= D(3)*st*beta - D(2 


*ct*(1+epsi) 


+ Hor; 


G(2,4) 


= 0; 






G(3,1) 


= G(1 ,3); 






G(3,2) 


= G(2,3); 






G(3,3) 


= D(2)*(1+epsi) A 2 + 


D(3)*beta A 2 


- Xi; 


G(3,4) 


= 0; 






G(4,1) 


= G(1 ,4); 






G(4,2) 


= G(2,4); 






G(4,3) 


= G(3,4); 






G(4,4) 


= D(1); 







%.... Form stiffness matrix and residual 
for i=1 :nbasis 

[hi,dhi] = basis(i,z,xlength,nbasis); mm = 3*(i-1); 
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for j=1 :nbasis 

[hj,dhj] = basis(j,z,xlength,nbasis); nn = 3*(j-1); 

% Compute noting the sparse structure of 

GB(1 :4,1 :2) = dhj*G(1 :4,1 :2); 

GB(1 :4,3) = hj*G(1 :4,3) + dhj*G(1 :4,4); 

BGB(1 :2,1 :3) = dhi*GB(1 :2,1 :3); 

BGB(3,1 :3) = hi*GB(3,1 :3) + dhi*GB(4,1 :3); 

% Assemble the result into the A matrix 

A(mm+1 :mm+3,nn+1 :nn+3) = A(mm+1 :mm+3,nn+1 :nn+3) + BGB*factor; 
end % loop on j 

% Form integral part of residual force and assemble into matrix 

b(mm+1:mm+3) = b(mm+1:mm+3) + (dhi*force + hi*c2)*factor; 

A(mm+1 :mm+3,ndm) = A(mm+1 :mm+3,ndm) - hi*c3*factor; 

end % loop on i 
z = z + dz; 
end % Loop on n 

%.. Add end load terms to the residual and coefficient matrix 
for i = 1 :nbasis 

[hi,dhi] = basis(i,1,xlength,nbasis); mm = 3*(i-1); 
b(mm+1:mm+3) = b(mm+1 :mm+3) - hi*c4*x(ndm); 

A(mm+1 :mm+3,ndm) = A(mm+1 :mm+3,ndm) - hi*c4; 
end 

%.. Add arc-length constraint terms to Hessian and residual 
b(ndm) = b(ndm) + dot(x-xo,x-xo) - alpha A 2; 

A(ndm,1 :ndm) = 2*(x-xo)'; 

%.. Transform Stiffness Matrix and Load vector terms to Global System 
A(1 :ndm-1 ,1 :ndm-1)=T'*A(1 :ndm-1 ,1 :ndm-1)*T ; 

A(1 :ndm-1 ,ndm)=T'*A(1 :ndm-1 ,ndm); 
b(1:ndm-1)=T'*b(1:ndm-1); 

return 



% 1 2 3 4 5 6 7 

% Evaluate ith basis function h and derivative dh BASIS 

function [h,dh] = basis(i,z,xlength,nbasis) 

L=1/xlength; 

%.. Quadratic Lagrangian finite element interpolation functions 
if nbasis==3 
if i==1 

h=1-3*z+2*z A 2; 
dh=(-3+4*z)*L; 
end 
if i==2 
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h=4*z-4*z A 2; 
dh=(4-8*z)*L; 
end 
if i==3 
h=-z+2*z A 2; 
dh=(-1+4*z)*L; 
end 
end 

%.. nbasis=5 
if nbasis==5 
if i==1 

h=1-25/3*z+70/3*z A 2-80/3*z A 3+32/3*z A 4; 
dh=(-25/3+140/3*z-80*z A 2+128/3*z A 3)*L; 
end 
if i==2 

h=16*z-208/3*z A 2+96*z A 3-128/3*z A 4; 
dh=(16-416/3*z+288*z A 2-512/3*z A 3)*L; 
end 
if i==3 

h=-12*z+76*z A 2-128*z A 3+64*z A 4; 
dh=(-12+152*z-384*z A 2+256*z A 3)*L; 
end 
if i==4 

h=1 6/3*z-1 1 2/3*z A 2-224/3*z A 3-128/3*z A 4; 
dh=(16/3-224/3*z+224*z A 2-512/3*z A 3)*L; 
end 
if i==5 

h=-z+22/3*z A 2-1 6*z A 3+32/3*z A 4; 
dh=(-1+44/3*z-48*z A 2+128/3*z A 3)*L; 
end 
end 

%.. nbasis=7 
if nbasis==7 
if i==1 

h=1 -1 47/1 0*z+406/5*z A 2-441 /2*z A 3+31 5*z A 4-1 1 34/5*z A 5+324/5*z A 6 ; 
dh=(-1 47/1 0+81 2/5*z-1 323/2*z A 2+1 260*z A 3-1 1 34*z A 4+1 944/5*z A 5)*L; 
end 
if i==2 

h=36*z-1566/5*z A 2+1044*z A 3-1674*z A 4+1296*z A 5-1944/5*z A 6; 
dh=(16-3132/5*z+3132*z A 2-6696*z A 3+6480*z A 4-11664/5*z A 5)*L; 
end 
if i==3 

h=-45*z+1053/2*z A 2-4149/2*z A 3+3699*z A 4-3078*z A 5+972*z A 6; 
dh=(-45+1053*z-12447/2*z A 2+14796*z A 3-15390*z A 4+5832*z A 5)*L; 
end 
if i==4 

h=40*z-508*z A 2+2232*z A 3-4356*z A 4+3888*z A 5-1296*z A 6; 
dh=(40-1016*z+6696*z A 2-17424*z A 3+19440*z A 4-7776*z A 5)*L; 
end 
if i==5 

h=-45/2*z+297*z A 2-2763/2*z A 3+2889*z A 4-2754*z A 5+972*z A 6; 
dh=(-45/2+594*z-8289/2*z A 2+11556*z A 3-13770*z A 4+5832*z A 5)*L; 
end 
if i==6 

h=36/5*z-486/5*z A 2+468*z A 3-1026*z A 4+5184/5*z A 5-1944/5*z A 6; 
dh=(36/5-972/5*z+1 404*z A 2-41 04*z A 3+51 84*z A 4-1 1 664/5*z A 5)*L ; 
end 
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if i==7 

h=-z+1 37/1 0*z A 2-1 35/2*z A 3+1 53*z A 4-1 62*z A 5+324/5*z A 6; 
dh=(-1 +1 37/5*z-405/2*z A 2+61 2*z A 3-81 0*z A 4+1 944/5*z A 5)*L; 
end 
end 

return 



% 1 2 3 4 5 6 7 

% Print and plot results of current step RESULTS 

function [] = 

results(x, indice, nelem, XX, YY,xlengthT,n,nu, test, ndm,nbasis,nnstep,npts,phiT) 
xp=zeros(3*nbasis,1); 
axis('square'); hold on; 
for nel=1 :nelem 
xlength=xlengthT(nel); 
phi=phiT(nel); 

%.. Transformation Matrix 
To=[cos(phi) sin(phi) 0; 

-sin(phi) cos(phi) 0; 

0 0 1 ]; 
for i=1 :nbasis 

T (3*i-2 :3*i ,3*i-2 :3*i)=To ; %complete T 

end 

for i=1 :3*nbasis 
if indice(nel,i) ~= 0 
xp(i)=x(indice(nel,i)); 
else 
xp(i)=0; 
end 
end 



%.. Transform x to local system 
xp=T*xp; 

%.. Determine if current step is an output step 
if mod(n,nnstep) == 0 

%.... Compute and print current geometry of beam 

z = 0; dz = 1/(npts-1); 
for ii = 1:npts 
u = zeros(3,1); 
for i = 1 :nbasis 

[h,dh] = basis(i,z,xlength,nbasis); mm = 3*(i-1); 
u(:> = u(:) + h*xp(mm+1 :mm+3); 
end 

yl(ii) = z*xlength + u(1); y2(ii) = u(2); 
z = z + dz; 
end 



for i=1:npts 
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Y1(i)=XX(nel)+y1(i)*cos(phi)-y2(i)*sin(phi); 

Y2(i)=YY(nel)+y1(i)*sin(phi)+y2(i)*cos(phi); 

end 

fprintf('%5i%9.5f%14.4e%14.4e%14.4e%4i%12.2e\n',n,x(ndm),u',nu,test) 

plot(Y1,Y2,'-'); 

end 

end 

return 
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RESUMEN 

Los materiales que se aplican para la remocion eficiente de arsenico (As) son los que 
tienen afinidad quimica con el metaloide y sus formas hidrolizadas. Las particulas de oxidos de 
hierro exhiben excelentes propiedades conducentes a la captura selectiva de contaminantes 
traza como el arsenico, pero por sus tamanos diminutos no pueden ser usadas en aplicaciones 
tipo flujo piston porque ocasionan elevadas caidas de presion. Encapsulando estas particulas 
dentro de materiales soporte mejora sus caracteristicas hidraulicas y mecanicas. 

En este proyecto se preparo un nuevo material sorbente mediante la aplicacion de un 
procedimiento quimico simple. El material sorbente hfbrido esta compuesto de una matriz de 
quitosano (biopolimero) y particulas de oxidos de hierro dispersadas dentro de las esferas del 
biopolfmero. El material una vez preparado fue usado para remover arsenico desde aguas 
contaminadas. Las pruebas de sorcion de arsenico usando agua sintetica muestran que lOpg/L 
de As se alcanzan luego de tratar 3000 volumenes de lecho de agua contaminada. Ademas, 
una planta prototipo cargada el material hfbrido puede tratar 12500 L de agua contaminada de 
origen natural sin exceder los 50 pg/L. El material puede tambien reusarse luego de su 
regeneracion con sosa caustica y acondicionamiento con agua acidificada. 

ABSTRACT 

Materials that are used for arsenic (As) removal are those with chemical affinity towards 
the metalloid and its hydrolyzed forms. Iron oxides particles exhibit excellent properties 
conducive to selective capture of trace contaminants such as arsenic. However, these tiny 
particles can not be used in plug flow applications due to the development of high pressure loss. 
Harnessing these particles within supporting materials improve their hydraulic and mechanical 
characteristics. 

In this study, we prepared a new sorbent material by a simple chemical procedure. The 
hybrid sorbent is formed by a matrix of chitosan (biopolymer) and iron oxide particles dispersed 
within of the biopolymer spheres. The sorbent once prepared was used to remove arsenic from 
contaminated waters. Arsenic sorption tests using synthesized water show that 10 pg/L of As is 
reached after treating 3000 bed volumes of contaminated water. In addition, a pilot plant loaded 
with the hybrid sorbent can treat 12500 L of natural water contaminated with arsenic without 
exceeding 50 pg/L. The material can also be reused after its regeneration with caustic soda and 
conditioning with acidified water. 



1 INTRODUCCION 



El arsenico (As) es un metaloide toxico que esta presente en las aguas superficiales y 
subterraneas en forma trivalente y pentavalente. La ingestion de agua contaminada con 
arsenico es un problema de salud publica a nivel mundial [1,2], pues por exposicion cronica, el 
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arsenico se acumula en el organismo y a ciertas concentraciones ocasiona afecciones en la piel. 
La evidencia epidemiologica de personas con ingestion prolongada de arsenico inorganico, via 
agua de bebida, es la hiperqueratosis palmo-plantar, cuya manifestacion principal es la 
pigmentacion de la piel y callosidades localizadas en las palmas de las manos y plantas de los 
pies. La presencia de arsenico en agua, su grado de contaminacion y la incidencia de 
enfermedades de piel en la India, Bangladesh, Taiwan, Vietnam, Argentina y Mexico se describen 
en diferentes estudios [3-11]. 

Por otra parte, las particulas de oxido de hierro son excelentes adsorbentes, de 
especies quimicas, incluyendo iones metalicos [12-14], aniones inorganicos [15,16] y algunos 
compuestos organicos [17,18]. Los oxidos de Fe(lll) hidratados (HFO) exhiben excelentes 
propiedades conducentes a la captura selectiva contaminantes traza por su comportamiento de 
acido debil. La predominancia de sus grupos funcionales y su capacidad de remocion depende 
del pH (pK a1 =6.5 y pK a2 =9) [19]. Por ejemplo a pH < 4 remueve arseniatos por un efecto 
combinado de interaccion electrostatica e interaccion acido-base de Lewis, en el rango de pH 
de 6.5 a 8.5 inmoviliza arsenico por interaccion acido-base de Lewis y a pH > 11 captura 
cationes metalicos [20]. Los metodos de preparacion de las particulas de oxidos de Fe(lll) son 
simples, sin embargo estas particulas diminutas no pueden ser utilizadas en columnas de lecho 
fijo, barreras reactivas, ni en ninguna configuracion tipo flujo piston, debido a que causan 
elevadas caidas de presion y porque no son durables. El encapsulamiento de las particulas de 
oxidos de Fe(lll) en materiales soporte mejora sus propiedades hidraulicas y mecanicas y las 
vuelve quimicamente mas estables. Investigaciones previas han utilizado como materiales de 
soporte arenas, alginato de sodio, carbon activado, etc., pero los resultados no han sido 
satisfactorios [20-24], 

En este estudio se prepararon esferas de quitosano, biopolfmero que se encuentra en 
forma natural en los caparazones de los crustaceos y en las paredes celulares de algunas 
plantas y hongos, que luego fueron impregnadas con particulas de oxido de Fe(lll) y utilizadas 
en la remocion de arsenico desde aguas sinteticas y naturales usando columnas de lecho fijo 
de laboratorio y en una planta piloto. 

2 METODOLOGIA 

2.1 Preparacion del material sorbente 

El procedimiento consta de las siguientes etapas: 

• Molienda y tamizado.- Granulos de quitosano son sometidos a molienda y tamizado para 
obtener tamanos de 250 pm. 

• Disolucion.- El quitosano tamizado es disuelto en una solucion de acido acetico al 4% y el 
producto se almacena a temperatura ambiente por 15 dias para disminuir la viscosidad. 

• Formacion de esferas.- La solucion se introduce en una jeringa de 50 mL y se deja caer gota 
por gota en una solucion al 5% de NaOH. 

• Lavado y entrecruzamiento.- Despues de 16 horas de contacto con la solucion alcalina, las 
esferas de quitosano son lavadas con agua desmineralizada y depositas en una solucion al 
2.5% de glutaraldehido. 

• Carga de oxidos de Fe (III).- 100 g de esferas de quitosano son sumergidas en un reactor 
conteniendo 500 mL de FeCI 3 al 3% y agitadas a 40 rpm por 48 horas. 

• Desorcion y precipitacion.- Las esferas extrafdas del reactor son sumergidas en 500 mL de 
una solucion de NaOH al 5% para la desorcion y precipitacion simultanea de los atomos de 
Fe(lll) en la superficie y poros de las esferas. 

• Lavado.- Las esferas impregnadas con los oxidos de hierro son sometidas a sucesivos 
lavados con agua desmineralizada para nivelar el pH a un valor igual a 7. 

2.2 Ensayos con columnas de lecho fijo 

Se condujeron usando una columna de 1.1 cm de diametro y 30 cm de longitud, 
empacada con quitosano impregnado con oxidos de Fe(lll) (Q-OFH). La razon diametro de 
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columna al diametro del sorbente fue de aproximadamente 10:1. Se corrieron dos series de 
columnas, la primera se alimento con agua sintetica preparada con 100 jug/L de arsenico, 120 
mg/L de sulfatos, 90 mg/L de cloruros y 100 mg/L de bicarbonatos y pH = 6.0. En la segunda 
serie se utilizo agua natural proveniente del rio Tuminguina. En las dos series de ensayos, se 
utilizo una bomba peristaltica de desplazamiento positivo, a un caudal de 1 mL/min. El efluente 
se recogio en tubos de ensayo colocados en un colector de fracciones Eldex PN 1243. 

2.3 Ensayos de regeneration 

En el proceso de regeneracion, se utilizo una solucion al 5% de NaOH. El regenerante 
se alimento a la columna de lecho fijo con un caudal de 0.5 mL/min. El efluente se recogio en el 
colector de fracciones y para posteriormente analisis. Las esferas y granulos regenerados 
fueron lavados con ocho volumenes de lecho de una solucion al 10% de acido acetico. 



2.4 Ensayos de digestion acida 

Para este ensayo, 250 mg de quitosano impregnado con OFH fueron sumergidos en 
dos lotes consecutivos de solucion acida (200 mL de H 2 S0 4 a 2N). Los lotes fueron colocados 
en un agitador rotativo a 40 rpm durante 48 horas. Los sobrenadantes fueron filtrados y 
analizados por hierro total. 

2.5 Analisis qui'micos 

Para la cuantificacion de arsenico en muestras tomadas del colector de fracciones, de 
la planta piloto y de las aguas del rfo Tuminguina, se utilizo el sistema Generador de Hidruros, 
Perkin Elmer, acoplado a un Espectrometro de Absorcion Atomica, Perkin Elmer, AA 100 (GH- 
AA100) equipado con una lampara de descarga de electrones a 197.3 nm de longitud de onda. 
Para medir el contenido de hierro total depositado en las esferas de quitosano se empleo 
igualmente el AA 100, Los cloruros, sulfatos y bicarbonatos fueron analizados siguiendo 
metodos estandares del Standard Methods. 

2.6 Reactivos qui'micos 

Todos los reactivos quimicos utilizados en esta investigacion fueron de grado analftico 
y comprados a la empresa Espectrocrom, representante en el Ecuador de Sigma-Aldrich 
Chemicals. El agua desmineralizada usada en este estudio, se la obtuvo desde un 
desmineralizador Millipore Direct-Q3. 

3 RESULTADOS Y DISCUSION 

3.1 Contenido de hierro total depositado en las esferas de quitosano 

El contenido de hierro total se cuantifico con la ayuda del espectrometro de absorcion 
atomica luego de una digestion acida usando el Metodo Fe-4500. Como se observa en la 
Figura 1 , el contenido de hierro en las esferas de quitosano varfa con el ciclo de carga. 

La Figura 1 indica que luego de cuatro ciclos se logra impregnar 38 mgFe/g de 
quitosano humedo. Notar que en los primeros ciclos, el contenido de hierro total es solamente 
22.3 mg/g, por lo que se requieren mas de dos ciclos de carga para tener una cantidad 
significativa de oxidos de hierro en el interior del quitosano. 

La carga de particulas de oxido de hierro en las esferas de quitosano es solamente 
moderada, debido que a pH acido (pH = 1.6) los grupos amino del quitosano se encuentran 
cargados positivamente y repelen a los atomos de Fe(lll). 
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Figura 1 Contenido de Fe total en las esferas quitosano en funcion del ciclo de carga 

3.2 Columna de lecho fijo usando agua sintetica 

En la Figura 2 se presenta el historial de la concentracion de arsenico en un ensayo de 
sorcion usando esferas de quitosano impregnadas con oxidos de Fe (III) y agua sintetica a pFH 
igual a 6.2. Notar que 50 pg/L se alcanzan luego de tratar 40,000 volumenes de agua 
contaminada con As. La selectividad por los oxianiones de As(V) puede explicarse porque a pH 
= 6.2, los sitios reactivos de los oxidos de Fe(lll) dispersos en las esferas de quitosano, 
predominante se encuentran cargados positivamente, por lo tanto remueven As mediante la 
combinacion de dos interacciones: i) electrostatica y ii) acido-base de Lewis y representadas 
por las reacciones [25]: 

I AR 

= FeOK> + H 2 AsO; ■ - - --> = FeOF£(H 2 AsO;) 

2(= FeOH^) + HAsO^" — l =^U(= FeOH£) 2 (HAsO^') 

LAB: interaccion acido base de Lewis 

Por otra parte, los grupos amino del quitosano a pH 6.2 son en su mayor parte 
positivos y en consecuencia retienen electrostaticamente a los oxianiones de As(V). Bajo 
condiciones ideales del proceso de sorcion, el porcentaje de grupos amino cargados 
positivamente responde a las siguientes ecuaciones [26]: 

QNH5->QNH 2 +H + 

Ka [QNH 2 ][H + ] 

[QNH 3 ] 

[QNH 2 ] t = [QNH5] + [QNH 2 ] 

De las ecuaciones 4 y 5, el porcentaje de grupos amino positivamente cargados puede 
obtenerse mediante la siguiente relacion: 

[QNH 3 ] _ 1 

= (6) 

[QNH 2 ] t 1 , Ka 

[H + ] 

Donde [QNH 2 ] representa el grupo amino neutro, [QNH 3 + ] es el grupo amino ionizado 
positivamente, [QNH 2 ] T constituye el grupo amino total presente en el quitosano y Ka es la 
constante de disociacion (mol/L). La pKa ha sido reportada por Roberts (1992b) en el rango de 

6.2 a 7.0 [27] 



(3) 

(4) 

(5) 



( 1 ) 

( 2 ) 
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Comparando con las arenas recubiertas con oxidos de hierro [21], el quitosano 
modificado muestra una mejor capacidad de sorcion porque los sitios reactivos estan no 
solamente localizados en la superficie externa sino tambien en el interior de la red de los poros 
de las esferas. 

Con los resultados de esta prueba se extrapolo la capacidad de remocion del quitosano 
(2734 pig As/g de quitosano humedo). Es preciso anotar sin embargo que el experimento se 
condujo con quitosano mixto (normal-desintegrado). De manera que los finos de las esferas 
pudieron generar una mayor area superficial y consecuentemente pudieron contribuir al 
incremento en la capacidad de remocion de As, aunque los finos no pueden utilizarse en filtros 
de plantas de tratamiento debido a que provocan excesiva cafda de presion. 
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Figura 2 Historial de concentraciones de As en una columna de lecho fijo 
empacada con quitosano modificado con oxidos de Fe(lll). 



Condiciones Experimentales 

Velocidad superficial = 0,6 m/h 
Tiempo de retencion = 1.52 min 



Afluente 

As (V) = 100 mg/L 
Cl' = 90 mg/L 
HC0 3 ' = 100 mg/L 
SO/' = 120 mg/L ^ i 
pH =6-2 4 



♦♦ 

♦ 



3.3 Columna de lecho fijo usando agua cruda del embalse del no Tuminguina 

Para evaluar la capacidad de remocion arsenico del quitosano impregnado con OFH se 
utilizo agua cruda del embalse de la EMAAP-Q, ubicado en el rfo Tuminguina. La concentracion 
de As T del agua reservorio fue de 49 gg/L y la del resto de componentes se indican en la Figura 
3. Antes del ensayo el agua natural fue ajustada a pH 6.0 y luego fue alimentada a la columna 
con una bomba peristaltica. Notar que la concentracion de 10 pgAsj/L, se alcanza luego de 
pasar 420 volumenes de lecho y que la concentracion del efluente de la columna es igual a la 
del influente (49 pg/L As T ) luego de pasar 8852 volumenes de lecho. En este ensayo ademas, 
se determino que la capacidad de remocion de As es igual a 1 79.7 gg As T /g sorbente, valor que 
comparado con el del ensayo con agua sintetica es de un orden de magnitud inferior. Este 
fenomeno puede ser atribuido a la competicion de la materia organica natural (40 mg/L, medida 
como carbono organico total) por los sitios de sorcion en los oxidos de Fe(lll) y a la formacion 
de complejos entre la materia organica natural y el As(lll), favoreciendo la elusion temprana del 
arsenico con el efluente de la columna. Redman et al. (2002), reportan que la presencia de 
materia organica provoca un desplazamiento y retraso de las especies arsenicales de los sitios 
de sorcion en sorbentes naturales, por la ocupacion u obstruccion de una gran proporcion de 
superficie reactiva, provocando de esta manera una disminucion en la cantidad de sitios activos 
para capturar arsenico[28]. 
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Figura 3 Historial de concentraciones de As en una columna de lecho fijo empacada 
con quitosano modificado con oxidos de Fe(lll) 



3.4 Columna de lecho fijo utilizando agua cruda con una pre-columna de carbon 
activado 

Con el proposito de evitar la rapida elusion de arsenico desde la columna de lecho fijo, 
se procedio a colocar una pre-columna de carbon activado para retener la materia organica 
natural. Los resultados de este experimento revelan que los 50 pg/L de As T , se alcanzan luego 
de pasar aproximadamente 10700 volumenes de lecho de agua natural contaminada con As y 
que se pueden remover 760 pg/L de As T (ver Figura 4). 
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Figura 4 Perfil de As en el efluente del arreglo experimental con dos columnas de lecho fijo 
usando agua de la vertiente de la laguna de Papallacta. 



Es evidente que la columna con carbon activado colocada antes de la columna con 
quitosano modificado produce un mejoramiento en la remocion de As desde el agua 
contaminada. Toles et al. (1997), senalan que el carbon activado remueve eficientemente 
compuestos organicos no polares y recalcitrantes en la potabilizacion de agua y en el 
tratamiento de aguas residuales [29], 

3.5 Regeneracion del quitosano impregnado con OFH 

Los materiales sorbentes son viables economicamente cuando estos son usados 
durante varios ciclos de sorcion y desorcion. Para examinar la factibilidad de que el Q-OFH sea 
reusado, se lo regenero con una solucion al 5% de NaOH. En la Figura 5, se muestra el perfil 
de concentraciones de arsenico total durante la regeneracion. Notar que con solo 10 
volumenes de sosa caustica se puede recuperar un 89% de arsenico retenido en el material 
durante la sorcion (673 pg As T /g quitosano). La observacion de que la desorcion de arsenico 
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fue completada en solo 10 volumenes de lecho, demuestra que los sitios de sorcion de las 
partfculas de OFH son facilmente accesibles. 
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Figura 5 Perfil de concentracion de arsenico durante la regeneracion del Q-OFH. 

A pFH alcalino, los sitios de sorcion en las partfculas OFH no contienen hidrogenos y se 
hallan cargados negativamente; De igual manera todas las especies arsenicales tienen carga 
negativa a este pH (aniones), de manera que el efecto de co-exclusion de Donnan es el que 
predomina dando como resultado una eficiente desorcion de las especies arsenicales [30]. 

La estequiometrfa de las reacciones de desorcion de los arseniatos y arsenitos 
adsorbidos en las superficies de las partfculas OFH, se expresan de la siguiente manera [25]: 



Regenerante: 

5% NaOH 




1 peO' 

/ 

O' 

X + 2AsO„ 3 ' + As0 2 ‘ + 1 1 H 2 0 



FeOhC 



✓FeOH 



H 2 As0 4 



HAs0 2 



paO' 

/ 

FeO' 



Complejos de esfera interna Exclusion de Donnan de Aniones 
pH < 8.0 



(7) 



Figura 6 llustracion del mecanismo de desorcion de las especies arsenicales 
durante la regeneracion del Q-OFH. 



Posterior a la regeneracion, las esferas de Q-OFH deben ser reacondicionadas con 0.3 
M de acido acetico para disminuir el pH. Se necesitaron solo ocho volumenes de solucion acida 
para que el pH alcance 5.5. Luego del reacondicionamiento, las esferas de Q-OFH estan listas 
para el siguiente ciclo de sorcion y ningun ajuste posterior de pH es necesario. En la Figura 7 
se resumen las reacciones qufmicas de protonizacion [25]. 
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pH >11,0 pH <8.0 



Figura 7 Procesos de acondicionamiento de Q-QFH 



3.6 Diseno, instalacion y operacion de una planta piloto para remover arsenico 
desde agua de la vertiente de la laguna de Papallacta 

Para el diseno de la planta se tomaron en cuenta los datos de los ensayos de sorcion 
del arsenico en el Q-OFH usando columnas de lecho fijo y en las pruebas cineticas de sorcion. 
La planta fue disenada para abastecer con agua segura a una comunidad de 20 personas, con 
un consumo promedio de de 15 litros agua/dfa/persona (agua para bebida y alimentacion). En 
los ensayos de sorcion se determino que el Q-OFH tiene una capacidad promedio de 0.5 
mgAs/g de quitosano humedo. Con este ultimo dato y los valores de velocidad superficial y 
tiempo de retencion, se dimensiono el filtro que contiene el material sorbente. Las 
caracteristicas del filtro se detallan en la Tabla 1. 



Tabla 1 Especificaciones del filtro de la planta piloto 



Especificaciones 

tecnicas 


Unidad 


Valor 


Caudal 


L/min 


0.50 


Altura 


cm 


45.00 


Diametro exterior 


cm 


21.60 


Diametro interior 


cm 


21.20 


Area transversal 


dm 2 


3.46 


Altura de lecho con Q-OFH 


cm 


14.44 


Densidad humeda del Q 


g/L 


1000.00 



Por otra parte, la dosificacion del cloruro ferrico usado en este estudio como 
coagulante, para inmovilizar el arsenico contenido en el regenerante usado, fue realizada en 
forma manual. En la Tabla 2 se enlistan los componentes de la planta piloto y el la Figura 8 se 
observa la distribucion de sus unidades. 



Tabla 2 Componentes de la planta piloto 



Componentes de la planta 
piloto 


Volumen 

L 


Tanque de almacenamiento de 


4000 


agua cruda 

Tanque para almacenar 


250 


regenerante 

Tanque para acondicionamiento 


250 


Tanque para almacenar solucion 


50 


estabilizadora de pH 
Tanque para recolectar los 


50 


precipitados de Fe(lll) 
Filtro de carbon activado 


10 



150 



En cambio en la Tabla 3 se indican las caracterfsticas de operacion de la planta. Las 
tres bombas de 0.25 HP que i) alimentan con agua cruda a la planta prototipo, ii) ayudan a 
circular la solucion regenerante por el filtro que contiene Q-OFH, iii) bombean la solucion acida 
luego de la regeneracion y acondicionan las esferas de Q/OFH a un pH igual a 6.0; fueron 
seleccionadas considerando un caudal de diseno (0.5 L/min). 




Figura 8 Planta prototipo para remover arsenico. 



Tabla 3 Caracterfsticas de operacion de la planta prototipo 


Parametro de Operacion 


Valor 


Contenido de As en la vertiente de la 
Laguna de Papallacta 
Masa As acumulada/dfa: c onsumo 


1 94 pg/L As 


diario de agua x concentracion de As 
Masa removido por quitosano: 
c apacidad quitosano x masa del 


58.2 mgAs/g/d 


quitosano 

Dfas antes de regeneracion: As 


2500 mg As 


removido/ Masa de As diaria 


42.96 


Tiempo para alcanzar 300 L/d 
Masa de Q-OFH: masa removido por 


10 h 


Q-OFH/capacidad del Q-OFH 
Volumen de Q-OFH: masa Q- 


5000 g 


OFH/densidad Q-OFH 


5 L 


Caudal de agua tratada: 
Velocidad superficial del agua: 


30 L/h 


caudal/area transversal 
Tiempo de retention: a Itura del Q- 


1 .44 cm/min 


OFH/velocidad superficial 


10 min 


Volumen regenerante (5% NaOH): 


50 L 



5.7 Primer ciclo de operacion de la planta prototipo 

En la Figura 9 se muestra la concentracion de As medido en el tanque de recoleccion 
del agua tratada durante el primer ciclo de operacion de la planta prototipo. 
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Figura 9 Concentracion de Ast en el primer ciclo de operation en la planta prototipo usando agua del 
embalse de la EMAAP-Q, localizado en el rfo Tuminguina. 



Aiju.nlel embdlse. lioTuminijuiii.i 

AsT = 194 ^g/L 
SO,* = 43.3 mg/L 
Cl - = 58.6 mg/L 
HC0 3 = 122 mg/L 
COT =49 mg/L 
Si0 2 =4.8 mg/L 
pH =6.8 






Conriiciones expeiimentales: 

Tiempo retention = 10 min 
Velocidad superficial =1.44 Umin 
Volumen empacado =5L 



Notar que luego de pasar 710 volumenes de lecho de agua cruda se alcanza la 
concentracion de 10 pg As T /L mientras que se pueden tratar mas de 2500 volumenes de lecho 
(12500 L) sin que se llegue todavia a los 50 pg/L de arsenico. La rapida elusion del arsenico 
desde el filtro conteniendo el Q-OFH puede estar asociada con la elevada concentracion del As 
en el agua cruda (194 pg/L). 

5.8 Primera regeneration del Q-OFH usado en la planta piloto 

En la Figura 10 se observa la regeneracion del Q-OFH empacado en el filtro de la 
planta prototipo. Fue necesario el uso de solo 10 volumenes de lecho (50 L de sosa caustica) 
para recuperar aproximadamente el 50% de arsenico desde el Q-OFH. 




Volumenes de lecho 



Figura 10 Perfil de concentracion de As durante la regeneracion del Q-OFH 
empacado en la planta piloto. 



4 CONCLUSIONES 

■ La impregnation con 38 mg Fe T /g quitosano humedo es bastante buena si se toma en 
consideration que a pH de la dosificacion con los atomos de Fe(lll) (pH = 1.6), los grupos 
amino del quitosano estan cargados positivamente y repelen los cationes de Fe. Notar que 
la dispersion de los OHF en una mayor area superficial favorece el incremento en la 
reactividad del quitosano y consecuentemente mejora la remocion selectiva de las especies 
arsenicales. 
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■ El quitosano impregnado con oxidos de Fe(lll) presenta varias ventajas con respecto a 
otros materiales soporte emergentes: i) disponibilidad de la materia prima para su 
procesamiento, ii) facil metodologia de preparacion, iii) buenas propiedades mecanicas e 
hidraulicas, iv) posibilidad de uso en varios ciclos de adsorcion y desorcion, v) no altera la 
calidad del agua porque no es toxico, requerimiento indispensable para ser utilizado en 
tratamientos de aguas de bebida, vi) eficiente regeneracion, requisito indispensable para 
que sea reusado y vii) biodegradacion rapida. 

■ El mecanismo de retencion ejercido por las esferas Q-OFH sobre las especies arsenicales, 
se fundamenta en la combinacion de interacciones electrostaticas y acido-base de Lewis. 
Los oxianiones arsenicales se ligan electrostaticamente con los grupos positivos de los 
oxidos de Fe(lll) mientras que la interaccion acido-base de Lewis forma complejos de 
esfera interna entre el arsenico (donador de un par de electrones) y los OFH (aceptador de 
par de electrones). 

■ El arreglo columna de pre-columna activado de carbon y columna de lecho fijo con Q-OFH 
exhibe una significativa mejora en la capacidad de remocion de As desde las aguas de la 
laguna de Papallacta. La pre-columna con carbon activado retiene de materia organica 
natural que compite por los sitios de sorcion con el arsenico, forma complejos con el As(lll) 
e impide la captura de arsenico. 

■ El volumen de agua tratada en la planta piloto (12500 L) en su primer ciclo de operacion es 
bastante halagador, considerando que la concentracion de As en el agua fue superior a 
190 pg/L. La eficiencia de tratamiento podrfa mejorarse si se modifica ligeramente el 
procedimiento de carga de los atomos de Fe(lll) en las esferas de quitosano y si se 
controla la salida y posterior deterioro del Q-OFH. 

■ La regeneracion del Q-OFH empacado en el filtro de la planta prototipo es moderada con el 
paso de solo 10 volumenes de sosa caustica. 

■ La dificultad de inmovilizar y almacenar las especies arsenicales disueltas en la sosa 
caustica luego de la regeneracion es el principal inconveniente que presentan las esferas 
de Q-OFH. 

RECOMENDACIONES 

■ Investigar otros procedimientos para la impregnacion de los oxidos de Fe(lll) en las esferas 
de quitosano a fin de obtener la misma o mayor cantidad de oxidos dentro del quitosano en 
menos ciclos de carga. 

■ Realizar mas pruebas con la planta prototipo para conocer la potencialidad de aplicarse en 
el tratamiento de agua contaminada con As en pequenos caserfos. 
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PETROQUIMICA EN EL ECUADOR 
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Centro de Investigaciones Cientfficas 
Escuela Superior Politecnica del Ejercito 

RESUMEN 



El desarrollo tecnologico en el area de hidrocarburos, es un reto que estan 
enfrentando los centros de investigacion superior. Con la finalidad de dar valor agregado a las 
reservas significativas de crudo. 

En el presente trabajo, se propone el diseno de un equipo piloto de destilacion 
fraccionada, formada por una caldera, columna empacada, unidad de reflujo, separador de 
muestras y condensador. 

La finalidad es separar en funcion del punto de ebullicion, los componentes de la 
gasolina comercial, para caracterizar los componentes por cromatografia de gases. 

Se inicia separando la fraccion de eter de petroleo, se analiza sus componentes 
hexanos, heptanos, octanos, y con esta informacion se propone aplicaciones para la 
industria nacional. 



ABSTRACT 



The technological development in the area of hydrocarbons is a challenge that has to 
assume the investigation centers. With the purpose of generating value added to the significant 
reservations of raw. 

Presently work intends the design and factory of a team pilot of fractional distillation, 
formed by a boiler, packed column, reflux unit, separator of samples and condensed. 

The purpose is to separate in function of the point of boil the components of the gasoline 
to characterize the components for Gas Chromatography. 

You proceed to separate the fraction of ether of petroleum and to analyze their 
components, hexanos, heptanos, and octanos and with this information to propose applications 
for the national industry. 

1 INTRODUCCION 

El petroleo es producto de los cambios naturales sobre la materia Organica, que se 
ha acumulado debajo de la superficie de la tierra en cantidades increfbles y ha sido 
descubierta y utilizada para satisfacer nuestras variadas necesidades de combustible. 

El Ecuador es poseedor de este recurso y se desea aportar con tecnicas que dan 
valor agregado a este recurso natural. 

Partiendo de hidrocarburos livianos como gasolina, diesel y empleando tecnicas de 
refinacion (destilacion) se ha logrado producir sustancias qufmicas que satisfacen 
necesidades industriales, sociales y de investigacion para sfntesis qufmica. 

Es importante refinar estos derivados para desarrollar tecnicas de separacion ffsico- 
qufmica, ademas proponer un diseno y construir con tecnologfa nacional el equipo de 
destilacion fraccionada, para caracterizar los componentes qufmicos de cada una de las 
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fracciones destiladas mediante la aplicacion de cromatografia de gas y proponer aplicaciones 
de estas fracciones para la industria. 

2 MATERIA PRIMA 

La gasolina es una mezcla de hidrocarburos derivada del petroleo que se emplea como 
combustible en motores de combustion interna con encendido a chispa. Tiene una densidad de 
720 g/L. Un litro de gasolina tiene una energta de 34,78 mega julios aproximadamente. Sin 
embargo, en terminos de masa, la gasolina tiene una energta de 48,31 MJ/Kg. . Se obtiene a 
partir de la nafta de destilacion directa, que es la fraccion Itquida mas ligera del petroleo 
(exceptuando los gases). La nafta tambien se obtiene a partir de la conversion de fracciones 
pesadas del petroleo (gasoil de vacto) en unidades de proceso denominadas FCC (craqueo 
catalttico fluidizado) o hidrocraqueo. La gasolina es una mezcla de cientos de hidrocarburos 
individuales desde C4 (butanos) hasta C11 como, por ejemplo, el metilnaftaleno. 

Normalmente se considera nafta a la fraccion del petroleo cuyo punto de ebullicion se 
encuentra aproximadamente entre 28 y 177 e C A su vez, este subproducto se subdivide en 
nafta ligera (hasta unos 100 S C) y nafta pesada (el resto). La nafta ligera es uno de los 
componentes de la gasolina, con numero de octano en torno a 70. La nafta pesada no tiene la 
calidad suficiente como para ser utilizada para ese fin, y su destino es la transformacion 
mediante reformado catalttica, proceso qutmico por el cual se obtiene tambien hidrogeno, a la 
vez que se aumenta el octanaje de dicha nafta. 

Ademas de la nafta reformada y la nafta ligera, otros componentes que se usan en la 
formulacion de una gasolina comercial son la nafta de FCC la nafta ligera isomerizada, la 
gasolina de pirolisis desbencenizada, butano, butanos, MTBE, ETBE, alquilato y etanol. Las 
formulas de cada refinerta suelen ser distintas (incluso perteneciendo a las mismas 
compahtas), en funcion de las unidades de proceso de que dispongan y segun sea verano o 
invierno. Si no esta refinada puede tener hasta 1000 ppm de azufre. Tiene alrededor de un 
40% de aromaticos y 20% de olefinas. Sus numeros de octano (MON/RON) estan en torno a 
80/93. 

3TEORIA DE DESTILACION 

En la mezcla simple de dos Itquidos solubles entre si, la volatilidad de cada uno es 
perturbada por la presencia del otro. En este caso, el punto de ebullicion de una mezcla al 
50%, por ejemplo, estarta a mitad de camino entre los puntos de ebullicion de las sustancias 
puras, y el grado de separacion producido por una destilacion individual dependerta solamente 
de la presion de vapor, o volatilidad de los componentes separados a esa temperatura. Esta 
sencilla relacion fue enunciada por vez primera por el qutmico trances Frangois Marie Raoult y 
se llama ley de Raoult. Solo se aplica a mezclas de Itquidos muy similares en su estructura 
qutmica, como el benceno y el tolueno. En la mayoria de los casos se producen amplias 
desviaciones de esta ley. Si un componente solo es ligeramente soluble en el otro, su 
volatilidad aumenta anormalmente. 

3.1 Destilacion fraccionada 

Si se consigue que una parte del destilado vuelva del condensador y gotee por una 
larga columna apilada, en la que el liquido fluye hacia abajo sobre una pila de anillos o trocitos 
de tubertas de vidrio se llama destilacion fraccionada. La unica desventaja de la destilacion 
fraccionada es que gran parte (mas o menos la mitad) del destilado condensado debe volver a 
la parte superior de la torre y eventualmente debe hervirse otra vez, con lo cual hay que 
suministrar mas calor, las columnas empacadas son usadas para destilacion. 



El contacto liquido - gas en una columna empacada es continua, no por etapas, como 
en una columna de platos. El flujo de liquido cae hacia abajo en la columna sobre el area de 
empaque y el gas o vapor asciende en contracorriente, en la columna. En algunas columnas de 
absorcion de gases se usa corrientes en flujo. El comportamiento de una columna empacada 
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depende mucho del mantenimiento de una buena distribution de Ifquido y gas a traves del 
lecho empacado y esto es una consideration importante en el diseno de columnas empacadas. 



Un diagrama esquematico se indica en la figura 1, en la cual se indican las partes de 
una columna empacada Una columna de destilacion con empaque sera similar a la columna 
de platos. 




Entra da del 




Entrada de Ifquido 



Platos de soporte 
Espaciador 



liquido 



Figura 1 Columna Empacada. 



3.2 Procedimiento de diseno de columnas empacadas 

El diseno de una columna empacada involucra las siguientes etapas: 

• Seleccionar el tipo y tamano de empaque. 

• Determinar la altura de la columna para la separacion especificada. 

• Determinar el diametro de la columna (capacidad), para manipular los flujos de Ifquido 
y vapor. 

• Seleccionar y disenar los dispositivos internos de la columna, soportes de empaque, 
distribuidos de liquido, redistribuidores. 



4 FUNDAMENTO TEORICO 

4.1 Tipos de empaque 

Los principales requerimientos de los empaques son: que estos deberfan proporcionar 
un area grande y una area interfacial entre el gas y el Ifquido tambien grande; Adicionalmente 
tener una estructura abierta, baja resistencia al flujo de gas , promover una distribucion 
uniforme de Ifquido sobre la superficie del empaque y Promover un flujo de gas o vapor 
uniforme a traves de la seccion transversal de la columna. 

Diversos tipos y formas de empaques han sido desarrollados para satisfacer estos 
requerimientos los cuales generalmente son divididos en tres clases. 
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4.1 .1 Empaques vaciados o al azar 

Estas son piezas discretas de empaques de una forma geometrica especifica las 
cuales son vaciadas o colocadas al azar en la columna y es el tipo de empaque que se aplica 
en el presente estudio. 

4.1 .2 Empaques estructurados o arreglados 

Estos se hacen en capas de malla de alambre u hojas corrugadas. Secciones o partes 
de estos empaques son colocados dentro de la columna. 

4.1.3 Parrillas 

Estas tambien son colocadas sistematicamente dentro de la columna. 

Los empaques al azar son las mas comunes en las practicas comerciales. Los 
empaques estructurados son menos comerciales, pero su mercado ha crecido rapidamente en 
la ultima decada. Las aplicaciones de las parrillas estan limitadas principalmente a servicios de 
transferencia de calor y lavado y/o donde se requieran altas resistencias a las incrustaciones. 

Aunque las columnas empacadas para absorcion y desorcion tambien pueden 
disenarse como un proceso por etapas, es usualmente mas conveniente usar la forma 
integrada de las ecuaciones diferenciales establecidas para las velocidades de transferencia de 
masa en un punto en la columna. 

Cuando la concentracion de soluto es pequena, (menor a 10 por ciento), los flujos de 
gas y liquido seran esencialmente constantes a lo largo de la columna y la altura de empaque 
requerida, Z, en terminos del coeficiente de transferencia total de la fase gas K a y la 
composicion del gas esta dada por: 

(D 

K 0 aP**y-y t 

Donde, G m es el flujo molar de gas por unidad de area de seccion transversal; K G es la 
composicion del gas; a es el area interfacial por unidad de volumen; p es la presion total; y 1 
y y 2 es la fraccion molar del soluto en el gas en el fondo y en el tope de la columna 
respectivamente; y e es la concentracion del gas en equilibrio con la concentracion del liquido en 
cualquier punto. 

En terminos del coeficiente de transferencia total de la fase liquida K L y la composicion 
del liquido: 

r*i dx 

Z= — (2) 

K z aC t x e - x 

Siendo, L m el Flujo molar de liquido por unidad de area de seccion transversal; a el area 
interfacial por unidad de volumen; C t es la concentracion molar total; x 1 y x 2 es la fraccion 
molar del soluto en el liquido en el fondo y en el tope de la columna, respectivamente; x e es la 
concentracion del liquido en equilibrio con la concentracion del gas en cualquier punto; y e es la 
concentracion del gas en equilibrio con la concentracion del liquido en cualquier punto. 

5 TEORIA SOBRE COLUMNAS DE SEPARACION 

La relacion entre la concentracion de equilibrio y la concentracion actual se muestra en la 
Figura 2. 
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Concentracion de soluto en el liquido 



Figura 2 Relaciones de concentracion en absorcion de gases. 

Para propositos de diseno, es conveniente escribir las ecuaciones (1 ) y ( 2 ) en 
terminos de las “Unidades de Transferencia” ( HTU)\ donde el valor de la integral es el numero 
de las unidades de transferencia y el grupo frente al signo de integral, el cual tiene las unidades 
de longitud, es la altura de una unidad de transferencia. 

Z = Hog N 0 g ( 3 ) 

Z = H ol N ol 



Donde H oa es la altura de una unidad de transferencia de la fase gas. 



Hog~ 



K G aP 



(4) 



Siendo N OG es el numero de unidades de transferencia de la fase gas 



N, 



r v ' dy 



OG 



* y - y e 



H 0 i_ es la altura de una unidad de transferencia de la fase liquida 



Hql = 



K,aC t 



(5) 



( 6 ) 



N ol es el numero de unidades de transferencia de la fase liquida 





(7) 



El numero de unidades de transferencia de la fase gas es frecuentemente mas 
conveniente expresado en terminos de la presion parcial del gas soluto. 



160 




dp 

P~Pe 



( 8 ) 



La relacion entre la altura total de una unidad de transferencia y el coeficiente individual 
de pelicula de las unidades individuales H L y H a , las cuales se basan en la fuerza impulsora 
dada por la concentracion a traves de las pelfculas de gas y de liquido, esta dada por: 

Hqg = Ha+m —r~ Hi 

L * (9) 

Hql = Hi + ~ Hq 

ML 



Donde m es la pendiente de la linea de equilibrio y GJL m la pendiente de la linea de 
operacion. 

El numero de unidades de transferencia se obtiene por integracion grafica o numerica 
de las ecuaciones ( 5 ), ( 7 ) y ( 8 ) 

Donde las Ifneas de operacion y equilibrio son rectas, y ellas pueden usualmente ser 
consideradas asf para sistemas diluidos. El numero de unidades de transferencia esta dado 
por: 



■Vi - 



Ay 
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Figura 3 Numero de unidades de transferencia Nog como una funcion 
de yi/y 2 con rnG m /L m como parametro. 



Donde A y, m es la media logaritmica de la fuerza impulsora (la concentracion), dada por: 

A yin = ^ ~ 

ln( / A y 2 ) 

A yi = yi- Ye * 
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A y 2 = y 2 - y e * 



Si la curva de equilibrio y las lineas de operacion pueden tomarse como lineas rectas y 
la alimentacion de solvente es esencialmente libre de soluto, el numero de unidades de 
transferencia es dado por: 



Nqg= 



1 

1 -(mGJLJ 



In 



/ 

1 

\ 



mG„ ) 

) 



y\ + mG M 
y 2 



Esta ecuacion es graficada en la figura 3 la cual puede ser usada para hacer un 
estimado rapido del numero de unidades de transferencia requerido para una separacion dada. 

De la figura. 3 se puede ver que el numero requerido de etapas para una separacion 
dada es muy dependiente de la velocidad de flujo L m . Si la velocidad de solvente no se fija por 
otras consideraciones del proceso, la figura 3 se puede usar para hacer estimados rapidos de 
la altura de la columna para diferentes velocidades de flujo para encontrar el valor mas 
economico. Coldburn (1939) sugirio que el valor optimo del termino mGn/L m debe estar entre 
0.7 y 0.8. 

5.1 Prediccion de la altura de una unidad de transferencia 

No hay un metodo enteramente satisfactorio para predecir la altura de una unidad de 
transferencia. En la practica el valor para un empaque particular dependera no solamente de 
las propiedades ffsicas y velocidades de flujo del gas y liquido, sino tambien de la uniformidad 
de la distribucion de liquido a traves de la columna, la cual depende de la altura de la columna 
y del diametro. Esto hace dificultoso extrapolar datos obtenidos a partir de pequenas columnas 
de laboratorio y de plantas piloto a tamanos de columnas industriales. Siempre que sea posible 
los estimados deben basarse en valores actuales obtenidos a partir de columnas en operacion 
de similar tamano a la que se esta disenando. Valores experimentales para varios sistemas 
son dados por Eckert (1963) y Connel et at. (1960). Valores para el termino K G a partir de los 
cuales se puede estimar la altura de la fase gas H OG , se pueden encontrar en diferentes 
publicaciones; Eckert (1975), Norman (1961), Perry y Green (1980). El termino K G a es 
normalmente usado cuando se reportan datos experimentales de los coeficientes de 
transferencia de masa para empaques, asi el area interfacial efectiva para la transferencia de 
masa sera menor que el area superficial actual del empaque. 



Muchas correlaciones se han publicado para predecir la altura de una unidad de 
transferencia, y los coeficientes de transferencia de masa. El metodo aplicado en este estudio 
es fiable para trabajos preliminares de diseno con un adecuado factor de seguridad. 



Tabla 1 Valores tlpicos de la HTU para empaques al azar 



Empaque 


HTU (UNIDADES DE 
TRANSFERENCIA) 


25 mm 


0,3 a 0,6 m 


38 mm 


0,5 a 0,75 m 


50 mm 


0,6 a 1 ,0 m 
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5.2 Experimental 

El primer paso efectuado en el laboratorio es la destilacion simple. En este proceso se 
aprovecha que cada compuesto tiene una temperatura caracteristica de ebullicion. En el 
laboratorio basta una manta calefactora, un matraz y un refrigerante para separar la mezcla 
por destilacion. Sin embargo, en la industria se requieren enormes equipos que se denominan 
torres de destilacion. 

El equipo de una destilacion simple se presenta en la figura 4 y de una torre de 
destilacion primaria se muestra en la figura 5, Para el diseno de una torre de destilacion es 
un reto clasico estimar <i,Que diametro proponer para la columna? ^Cuantos platos de contacto 
Ifquido-vapor hacen falta o emplear rellenos y que tipo de relleno? ^En que lugar de la 
columna se alimenta con el hidrocarburo? que temperatura debe introducirse? ^Cuanto 
vapor alimentar al rehervidor? 




Figura 4 Destilacion simple de gasolina 
Tabla 2 Temperatura de ebullicion de hidrocarburos 



Nombre 


Intervalo de 
temperatura de 
ebullicion (C) 


Numero de 
carbonos 


Uso 


Eter de petroleo 


20-80 


5 a 7 


disolvente 


Gasolina 


35 - 220 


5 a 12 


combustible de autos 


Querosin 


200 -315 


12a 16 


combustible de 
aviones 



La temperatura cambia a lo largo de la columna. En la parte superior se tiene la mas 
baja temperatura, donde se encuentran en equilibrio los componentes mas ligeros (y de menor 
punto de ebullicion). Por el contrario, en la parte inferior la temperatura es mucho mas alta y es 
tambien la proporcion de los componentes pesados y menos volatiles. 

Como se colectan diversas fracciones de la columna de destilacion, la gasolina logra 
separarse en varias componente, cada una con un diferente intervalo de temperaturas de 
ebullicion e hidrocarburos de diferente numero de carbonos en su cadena (vease la tabla 2). 
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Figura 5 Columna de destilacion fraccionada experimental. 



La mayor parte de la gasolina que se destilan directamente de los crudos son n 
parafinas en forma predominante, estos materiales tienen propiedades antidetonantes como 
ejemplo tenemos al Hexano alrededor de 2%, n Heptano. 

La volatilidad se estudia de acuerdo a la curva de destilacion ASTM. La volatilidad de 
una gasolina se defina como la tendencia a pasar a fase vapor en unas condiciones 
determinadas. Para una gasolina concreta nos interesa conocer el IBP (PIE): punto inicial de 
ebullicion y el PFE: punto final de ebullicion, punto en el que se recoge: 

■ 10%destilado 

■ 20%destilado 

■ 50%destilado 

■ % total de residuo 

■ % total perdidas 



En el BOE, para la gasolina 97 nos interesa, ademas, conocer los siguientes datos: 

■ % recodido a 70 a C 

■ % recogido a 100 a C 

■ % recogido a 180 e C 



Como una regia general para empaques de ceramica puede tomarse como el 25 por 
ciento del peso del empaque. 
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Experimentalmente se desarrollan someramente las caracterfsticas de la cromatograffa 
gaseosa, enfocada hacia el analisis de las mezclas de hidrocarburos. Por otra parte, se 
menciona que el analisis cromatografico de una mezcla gaseosa proporciona los datos 
necesarios para calcular, con razonable exactitud, casi todos los parametros de interes de 
dicha mezcla. Esto es asi gracias a dos factores concurrentes: 

• El cromatograma de una mezcla de gases proporciona la composicion individual 
de todos los componentes mayoritarios de la mezcla. 

• Las propiedades de los gases son basicamente aditivas. En otras palabras, las 
propiedades de la mezcla son calculables a partir de las propiedades de los 
componentes individuales. La composicion compleja de los componentes 
puede apreciarse en la figura 6. 






Figura 6 Cromatograma de hidrocarburos 



El cromatograma fue obtenido mediante una columna de alta resolucion (columna 
capilar) y un detector FID, que solo es sensible a los hidrocarburos de la mezcla. 

• La serie de parafinas normales (los picos mas altos y casi equi-espaciados en el 
cromatograma, es perfectamente identificable. 

• Existe un “fondo” importante de componentes no resueltos (separados en forma 
incompleta) correspondientes a compuestos ramificados, aromaticos, cfclicos, 
etc. Sin embargo, pese a la complejidad aparente de este cromatograma, el 
reporte de laboratorio que acompana los estudios PVT suele adoptar la siguiente 
forma. 



Tabla 3 Porcentaje molar 


Componente 


[% Molar] 


i-Pentano 


1.463 


n-Pentano 


2.697 


Hexanos 


5.037 


Heptanos 


7.569 
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Octanos 


8.689 


Nonanos 


7.171 


Decanos 


6.589 


Undecanos 


5.327 


Dodecanos 


4.359 


Tridecanos 


4.403 


Tetradecanos 


3.621 


Pentadecanos 


2.830 


Hexadecanos 


2.472 


Heptadecanos 


2.434 


Octadecanos 


2.121 


Nonadecanos 


1.903 


Eicosanos y sup. 


27.216 



Esta informacion, que aparece en la tabla 3 que indica los porcentajes molares de 
hidrocarburos, es la que se emplea como datos de entrada en los estudios de simulacion de 
termodinamica. 



6 CONCLUSIONES 

• En el presente articulo se presenta el diseno de una columna de fraccionamiento a la 
escala piloto. 

• Se estima que la composicion de la muestraquees una mezcla liquida de 
hidrocarburos no se mide en forma directa, sino que es el resultado de una serie de 
mediciones, calculos y simplificaciones que conducen a los valores. 

• La caracterizacion de las diferentes fracciones “livianas” se obtendra mediante 
cromatograffa de alta resolucion. 

• Las propiedades de la fraccion “pesada” se obtendra por una via independiente 
(preferentemente mediante una destilacion fraccionada. 

• La altura de la columna de 80 cm. se calculo para que la relacion optimo se encuentre 
entre mG„/Z. m = 0,6 a 0,8. 

• Como una gufa general, la altura maxima del lecho no debera exceder a 3 tantos el 
diametro de la columna para anillos Rasching, y 8 a 10 para anillos Pall y monturas 
Intalox y se esta ensayando para empaques vaciados al azar entre 13cm y 40 cm. 

Para un diametro de la columna de 4.5 cm. 

• Es de notar que existen numerosas herramientas de control que permiten establecer la 
contiabilidad de la composicion obtenida como densidad, tndice de refraccion, punto de 
ebullicion, viscosidad. 

• Es necesario un estimado de la cantidad de liquido retenido en el empaque y esta fue 
ensayada experimentalmente cuya resultado de retencion es de 15 ml para una altura 
de 40 cm. de empaque. 
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• La destilacion se basa en la transferencia de masa entre las fases liquido - vapor de 
una mezcla de hidrocarburos. Cuya relacion se monitorea mediante la medicion del 
flujo de la materia prima y del destilado. 

• La destilacion permite la separacion de los componentes de una mezcla de 
hidrocarburos, como es la gasolina, en funcion de sus temperaturas de ebullicion. 

• Para que se produzca la “separacion o fraccionamiento” de los codes, se debe alcanzar 
el equilibrio entre las fases liquido-vapor, ya que de esta manera los componentes mas 
livianos o de menor peso molecular se concentran en la fase vapor y por el contrario los 
de mayor peso molecular predominan en la fase liquida, en definitiva se aprovecha las 
diferencias de volatilidad de los hidrocarburos. 

• El equilibrio liquido-vapor, depende principalmente de los parametros termodinamicos, 
presion y temperatura del sistema. Las unidades se disenan para que se produzcan 
estos equilibrios en forma controlada y durante el tiempo necesario para obtener las 
fracciones especificadas. 

• Basicamente el proceso consiste en vaporizar los hidrocarburos de la gasolina y luego 
condensarlos en codes definidos. Modificando fundamentalmente la temperatura, a lo 
largo de la columna fraccionadora. 

• La fase liquida se logra con reflujos o reciclo de hidrocarburos retornados a la torre. 
Estos reflujos son corrientes liquidas de hidrocarburos que se enfrfan por intercambio 
con crudo o fluidos refrigerantes. La funcion u objetivo principal de estos, es eliminar o 
disipar en forma controlada la energfa cedida a los hidrocarburos en la caldera, de esta 
manera se enfrfa y condensa la carga vaporizada, en codes o fracciones de 
hidrocarburos especfficas, obteniendose las fracciones correspondientes. 

• En la presente investigacion se emplea como materia prima gasolina extra de libre 
venta. 

• Un equivalente barato de la torre de burbujeo es la llamada columna apilada, en la que 
el liquido fluye hacia abajo sobre una pila de empaque de vidrio. 

• La unica desventaja de la destilacion fraccionada es que gran parte (mas o menos la 
mitad) del destilado condensado debe volver a la parte superior de la torre y 
eventualmente debe hervirse otra vez, con lo cual hay que suministrar mas calor. Por 
otra parte, el funcionamiento continuo permite grandes ahorros de calor, ya que el 
destilado que sale puede ser utilizado para precalentar el material que entra. 
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CARACTERIZACION ESPACIAL DE HABITATS EN UN 
PAISAJE FRAGMENTADO DE SANTO DOMINGO 
DE LOS TSACHILA 

Jaime Villaci's, Diego Cortez, Jorge Zurita, Marco Luna, Carlos Chiriboga. 

Departamento de Ciencias de la Vida 
Escuela Superior Politecnica del Ejercito 

RESUMEN 

El objetivo general de este estudio fue generar informacion que ayude a enriquecer el 
conocimiento sobre la estructura y composicion de paisajes fragmentados y conocer como los 
ecosistemas forestales contribuyen a la conservacion de la diversidad arborea de la zona. Para 
esto se utilizo una fotografia aerea, la cual se la digitalizo a escala de 1 :5000 en ambiente SIG 
y se utilizo la guta de fotointerpretacion propuesta por el proyecto FRAGMENT. 

La composicion y la configuracion espacial de los habitats en el paisaje, se determino 
cuantificando la cantidad de parches correspondiente a cada clase o categorta de habitat. Para 
la caracterizacion arborea de los habitats en el paisaje, se seleccionaron dos habitats (bosque 
secundario y bosque ripario) y se establecieron 8 parcelas de muestreo por cada habitat. 

La informacion se analizo con estadtstica descriptiva, ANDEVAS y pruebas de 
comparacion de medias de Duncan (5%). Se determino el uso del suelo en 23275,4497 ha y 
se identificaron 9 habitats diferentes. El habitat potrero con arboles ocupo el 40,70% del area 
de estudio, seguido de cultivos de ciclo corto y cultivos perennes que alcanzaron el 14,58% y 
25,20% respectivamente. Se encontraron 60 especies de arboles en los dos habitats 
inventariados (en el bosque ripario 40 especies y en el bosque secundario 43 especies). 

Se encontraron 27 familias botanicas diferentes; la familia Fabaceae - Mim fue la mas 
representativa ya que agrupo al 17,32% de las especies Inventariadas, mientras que el 11 ,29% 
del total de especies inventariadas se agruparon en 15 familias. El 16,19% de los arboles 
inventariados corresponden a “Pechuga de gallina” (Aegiphila alba), seguido por “Pachaco” 
( Schizolobium parahyba) que agrupo al 9,22% del total de especies. 

De acuerdo a los resultados del Indice de Valor de Importancia (IVI%) se observo que 
las especies con mayor valor fueron: Aegiphila alba con 12,062% y Pouteria capacifolia con 
11,342%. El bosque secundario presento similar abundancia, riqueza y diversidad que el 
bosque ripario, lo que implica que se pueden establecer estrategias de conservacion similares 
en estos dos tipos de habitats para contrarrestar los efectos de la fragmentacion del paisaje 
sobre la biodiversidad. 
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1 INTRODUCCION 

En paisajes fragmentados es comun encontrar remanentes de bosques primarios inmersos en 
una matriz dominada por areas agrfcolas y ganaderas, donde se pueden encontrar diferentes tipos de 
cobertura arborea tales como bosques riparios, vegetacion secundaria de distintas fases de sucesion, 
arboles aislados, y arboles en Ifnea (Kleinn 2000). Esta cobertura arborea puede ser de suma 
importancia para la conservacion de la biodiversidad, en primera instancia por que forma parte de ella y 
ademas por que puede proveer de habitats y recursos indispensables para la sobrevivencia de los 
organismos en areas abiertas (Galindo- Gonzalez et al. 2000, Morales y Klein 2001). 

En anos recientes ha aumentado el interes por conocer la biodiversidad en areas influenciadas 
por actividades agrfcolas y ganaderas, con el fin de recomendar actividades de manejo para la 
conservacion (Villacfs 2003). Por otro lado, conocer la riqueza diversidad y composicion de arboles en 
areas abiertas, es vital para la conservacion de biodiversidad puesto que se ha demostrado su rol para la 
sobrevivencia de especies en lugares donde la intervencion humana es frecuente (Guevara et al 1998). 
Poco se conoce aun sobre la distribucion, composicion y estructura de los diferentes habitats que 
conforman los paisajes fragmentados tropicales y sobre la diversidad arborea que se puede encontrar en 
estos lugares. Conocer la composicion y estructura del paisaje es importante para la conservacion pues 
es una herramienta que permite entender a escalas mayores los distintos problemas ambientales 
relacionados con los organismos y documentar dicha informacion para estudios posteriores (Guariguata 
y Kattan 2002). 

La zona de Santo Domingo de los Tsachilas esta dotada de una alta biodiversidad, en una 
variedad de ecosistemas que se extienden en una zona con un alto fndice de especies endemicas, 
considerada por cientfficos ambientales como laboratorio para la investigacion genetica mundial. Por su 
ubicacion geografica, tiene un comercio muy rico ya que se ha convertido en un puerto terrestre de 
intercambio entre Sierra y Costa. La actividad mas importante es la ganaderfa, pues es aquf en donde se 
desarrolla el mayor mercado ganadero del pafs. El proposito de este estudio fue generar informacion que 
ayude a enriquecer el conocimiento sobre la estructura y composicion de paisajes fragmentados y 
conocer como los ecosistemas forestales contribuyen a la conservacion de la diversidad arborea de la 
zona de Santo Domingo de los Tsachilas. 

2 METODOS 

El presente estudio se realizo en un area de 23275 ha. El area de estudio se encuentra en la 
provincia de Santo Domingo de los Tsachilas, Region Litoral, latitud 0024‘5‘ S y longitud 79°I6‘36‘“ W. 
La zona se ubica a 300 m de altitud con una precipitacion media anual de 3045,1 mm. La humedad 
relativa promedio es del 85.8 % y la temperatura promedio es de 24,2 e C con poca variacion durante el 
dfa (INIAP 2004). 

Se caracterizo el paisaje en terminos de la composicion y estructura (area y distribucion) de los 
diferentes habitats (o usos del suelo) presentes en la zona. La fotograffa aerea se digitalizo a escala de 
1 :5000 en ambiente SIG (archivo tipo .dwg) y se utilizo la “gufa de fotointerpretacion” propuesta por el 
proyecto FRAGMENT (Morales y Kleinn 2001) para interpretar la fotograffa y determinar los habitats 
presentes en la zona. Los habitats se determinaron recorriendo el area de estudio, donde se observo 
diferentes sistemas de produccion pecuaria y cultivos, por lo cual se definio a consideracion personal los 
siguientes: cultivos de ciclo corto, cultivos perennes, potreros con arboles, potreros sin arboles, 
plantaciones forestales, bosques riparios, bosques secundarios y zonas pobladas. La definicion de los 
habitats mencionados, agruparon a todos los diferentes cultivos y formas de produccion ganadera de la 
zona. En cada habitat se tomo puntos de georeferenciacion con los GPS, con los cuales se identified con 
exactitud en la fotograffa aerea los colores que identificaron a cada habitat. Las medidas utilizadas en el 
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analisis del paisaje calculados (Kramer 1997, Coppedge et al. 2001 , Gallego 2002) fueron las siguientes: 
area total de la clase (ha), proporcion en el paisaje (%), numero de parches y densidad de parches 
(#/1 00 ha). 

Para analizar la composicion y la configuracion espacial (estructura) de los elementos o habitats 
en el paisaje, se determine la cantidad de parches correspondiente a cada clase o categorfa de habitat. 
Se estimo el area y la proporcion de cada uno de los tipos de habitats, ademas del numero de parches 
segun habitat y su densidad de parches en el paisaje. 

2.1 Composicion, estructura y diversidad arborea de diferentes habitats. 

2.1.1 Muestreo 

De los datos suministrados de la fotointerpretacion, se elaboro un mapa de uso del suelo o de 
habitats que se comprobo en el campo y fue la base para dirigir el proceso de muestreo y 
establecimiento de las parcelas para la caracterizacion arborea de los habitats en el paisaje. En total se 
ubicaron y establecieron 1 6 parcelas de muestreo (8 por cada habitat) de 1 00 m x 1 0 m y 20 x 50 m (0.1 
ha) en 2 de los habitats encontrados donde se visualizo un mayor numero de arboles (por razones 
economicas). Para ubicar estas parcelas, se seleccionaron al azar 16 puntos sobre el mapa de habitats 
(8 puntos por tipo de habitat) y mediante un sistema de posicionamiento global (GPS Garmin) se localizo 
cada uno de estos puntos en el campo. Este punto se denomino esquina principal y de allf se partio para 
establecer y orientar la parcela, en sentido norte-sur. En cada una de estas parcelas de 0,1 ha se 
midieron todos los arboles con un diametro a la altura del pecho mayor o igual a diez centfmetros (dap > 
10 cm), el DAP, la altura total el arae basal y el volumen total. Todas las especies de arboles fueron 
identificadas en el campo con la ayuda del Dr. Fernando Hurtado, especialista en botanica taxonomica. 

2.1 .2 Riqueza y diversidad arborea 

Para cada parcela se calculo el numero de individuos, y la riqueza (S), (que se refiere al numero 
de especies total presentes en un area determinada). La diversidad florfstica se analizo mediante el 
calculo de los indices de Shannon, Simpson (Magurran 1988), Berger - Parker, Margaleff, Alpha y 
Mackintosh. 

2.1 .3 Composicion y estructura arborea 

Se determino la abundancia, la frecuencia de individuos por especie y frecuencia de especies 
por familia para cada una de las 8 parcelas de cada habitat. Se calculo tambien el fndice de valor de 
importancia (IVI), para determinar las especies mas importantes por parcela y para cada uno de los 
diferentes habitats analizados. Para comparar la abundancia, diversidad, riqueza y dominancia entre 
bosque secundario y bosque ripario se realizaron ANDEVAS y pruebas de comparacion de medias de 
Duncan al 5 % (Infostat 2007). 

3 RESULTADOS Y DISCUSION 
3.1 Uso del suelo 

Se determino 7 diferentes habitats en la zona de estudio: cultivos de ciclo corto, potreros sin 
arboles, potreros con arboles, cultivos perennes, bosques, bosques riparios y plantaciones forestales 
(Cuadro 1). 
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Cuadro 1 Description de los diferentes tipos de habitats identificados en Santo 
Domingo de los Tsachilas Ecuador, 2008. 



2 HABITAT 


DESCRIPCION 


Cultivos de 
ciclo Corto 


Areas caracterizadas por presencia de especies de ciclo de vida inferior a un 
ano. Facilmente diferenciadas por colores claros mas intensos en 
comparacion con las areas de potreros. 


Potreros 


Areas dedicadas exclusivamente a la ganaderfa, identificados por presentar 
colores claros. 


Potreros con 
arboles 


Areas dominadas por arboles agrupados en las pasturas. En el area de 
estudio se identified gran cantidad de estas areas. Muestra colores claros al 
igual que en areas de potreros con la diferencia de pequenas manchas de 
colores oscuros dispersos en el area. 


Cultivo 

perenne 


Areas caracterizas por presencia de especies de ciclo largo, es decir superior 
a los dos anos. Identificados por colores oscuros y en la mayorfa de los 
casos alineados sistematicamente. 


Bosque 


Areas donde los arboles son la especie vegetal dominante. Se caracterizaron 
por la presencia de arboles maduros de diferentes edades, especies y porte 
variado. Presenta colores oscuros intensos sin ningun tipo de alineacion. 


Bosque 

Ripario 


Franjas alargadas de bosque secundario, encontradas en arroyos, rfos y 
esteros. 


Plantacion 

Forestal 


Areas de cultivos lenosos especialmente maderables y facilmente 
distinguibles debido a que se mostraban alineados sistematicamente. 



El area total de la zona de estudio fue 23275,4497 ha. El 14,57% del area total fue destinada a la 
siembra de cultivos de ciclo corto (mafz, yuca, malanga), el 25,19% a la siembra de cultivos perennes 
(cacao, citricos, palmito, palma africana, banano), el 40,70% a potreros con presencia de arboles, el 
10,42% a potreros sin arboles y el 0,07% a plantaciones forestales. En cuanto a bosques primarios y 
secundarios se destino un 7,16% del area total y a bosques riparios se destino el 0,96% (Cuadro 2). 

Cuadro 2 Representation porcentual de cada uno de los habitats en el paisaje fragmentado de Santo Domingo de 
los Tsachilas, Ecuador 2008. 



Habitat 


Total ha 


% del 

paisaje 


N°de 

parches 


Bosque 


1668,65 


7,17 


79 


Bosque ripario 


224,99 


0,97 


24 


Cultivos de ciclo corto 


3392,82 


14,58 


1632 


Potreros con arboles 


9474,08 


40,70 


2331 


Potreros sin arboles 


2426,90 


10,43 


436 


Cultivos Perennes 


5865,34 


25,20 


847 


Plantacion Forestal 


16,46 


0,07 


5 


Indeterminado 


111,20 


0,48 


8 


Poblados 


105,32 


0,45 


3 



El habitat que predomino en la zona fue potrero con arboles con 9474, 07 ha, seguido por 
cultivos perennes con 5865,33 ha, esto se debe a que en la mayor parte del area de estudio la topograffa 
de la zona dificulta la explotacion agrfcola, por lo que son tierras aptas para la siembra de pastos y 
produccion ganadera. El habitat potrero con arboles estuvo formado por 2331 parches, seguido por 
cultivos de ciclo corto que conto con 1632 parches, los habitats bosque ripario y bosque secundario 
estuvieron formados por 24 y 79 parches respectivamente. 
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Estos resultados son diferentes a loa encontrados por Chacon (2003) en Costa Rica, en donde 
encontro 12 tipos de habitats en 293 parches: bosque denso, vegetacion secundaria joven, bosque 
ripario, palmito, grupo de arboles en pasturas, plantaciones forestales y frutales, huertos, pastos y suelo 
desnudo. Tambien difieren de otro estudio realizado en Guanacaste, Costa Rica, en donde se 
determinaron seis tipos de habitats que dominaron el paisaje: fragmentos de bosque seco, bosques 
riparios, charrales, cercas vivas y potreros de alta cobertura de arboles (Cardenas 2002). 

3.2 Caracterizacion del componente arboreo 

Se encontro 60 especies de arboles en los dos habitats inventariados (en el bosque ripario 40 
especies y en el bosque secundario 43 especies). El numero total de arboles fue de 531 agrupados en 
27 familias botanicas. La altura promedio de los arboles fue de 14.85+0.45 m, el area basal promedio fue 
de 0.08+0.01 m 2 y el volumen promedio fue de 0.71 ±0.08 m3. 

Se agruparon 348 arboles de 10 especies, el 16,19% de los arboles inventariados 
correspondieron a “Pechuga de gallina” ( Aegiphila alba), siendo la especie de mayor importancia. Para la 
especie “Pechuga de gallina” se inventario 86 individuos, seguido por “Pachaco” con 49 individuos. El 
9,22% del total de individuos correspond^ a “Pachaco” ( Schizolobium parahyba) (Cuadro 3). Se 
encontro 27 familias botanicas diferentes en el inventario realizado. La familia Fabaceae - Mim fue la 
mas representativa ya que agrupo al 17,32% de las especies inventariadas. El 11,29% del total de 
especies inventariadas se agruparon en 15 familias. 



Cuadro 3 Especies mas comunes, nombre vulgar, familia, porcentaje y 
numero de parcelas con la especie. 



Nombre vulgar 


Especie 


Familia 


# 

individuos 


% 


# 

parcelas 
con la 

especie 


Pechuga de 

gallina 


Aegiphila alba 


VERBENACEAE 


86 


16,20 


12 


Pachaco 


Schizolobium 

parahyba 


FABACEAE 


49 


9,23 


4 


Chisparo 


Pithecellobium 

exelsum. 


FABACEAE 


48 


9,04 


4 


Inga 


Inga marginata 


FABACEAE 


36 


6,78 


5 


Caucho 


Castilla elastica 


MORACEAE 


33 


6,21 


12 


Visola 


Wettinia quinaria 


ARECACEAE 


33 


6,21 


6 


Zapotillo 


Pouteria capacifolia 


SAPOTACEAE 


20 


3,77 


11 


Nacedera 


Erythrina smithiana 


FABACEAE 


17 


3,20 


2 


Uva 


Pouroma bicolor 


CECROPIACEAE 


15 


2,82 


7 


Laurel 


Cordia alliodora 


BORAGINACEAE 


11 


2,07 


4 



De acuerdo a los resultados del indice de Valor de Importancia (IVI%) se observo que las 
especies con mayor valor fueron: Aegiphila alba (“Pechuga de gallina”) con 12,062% y Pouteria 
capacifolia (“Mamey”) con 11,342%. Las dos especies se encontraron distribuidas ampliamente en 
bosque ripario y bosque secundario. El 66,028% del IVI% correspond^ a 10 especies y el 33,97% a 51 
especies. El menor porcentaje de IVI% pertenecio a las especies Brownea multijuga (“Clavelin”) y Grias 
multinervia (“Grias”) con un IVI% de 0,1 1 8% (Cuadro 4). 



Cuadro 4 Abundancia relativa de la especie (A%), dominancia relativa de la especie (D%), frecuencia relativa de la 
especie e indice de valor de importancia (IVI%) de especies inventariadas mas comunes. 
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Nombre vulgar 


A% 


D% 


F% 


IVI % 


Pechuga de gallina 


16,20 


7,99 


0,89 


25,09 


Pachaco 


9,23 


6,09 


0,28 


15,59 


Chfsparo 


9,04 


10,11 


0,27 


19,42 


Inga 


6,78 


7,69 


0,34 


14,81 


Caucho 


6,21 


5,11 


0,80 


12,13 


Vfsola 


6,21 


5,11 


0,40 


11,72 


Mamey 


3,77 


18,98 


0,71 


23,46 


Nacedero 


3,20 


1,23 


0,13 


4,56 


Uva 


2,82 


2,33 


0,45 


5,61 


Laurel 


2,07 


2,00 


0,26 


4,33 



3.3 Distribucion de diametros a la altura de pecho y alturas en arboles inventariados. 

En las parcelas se encontraron un total de 519 arboles con dap > 10 cm. El 44,63% de los 
arboles inventariados presentaron un dap > 10 cm y < 20 cm. Tres arboles presentaron un dap > 90 cm 
que corresponde al 0,56%. El 96,40% de los arboles inventariados tienen dap < 50 cm que corresponde 
a 496 arboles (Figura 1). 




Figura 1 Distribucion de clases diametricas de arboles presentes en parcelas (n = 531 individuos), Santo Domingo 

de los Colorados, Ecuador, 2007. 



Se determino que el 29,75% de los arboles inventariados alcanzaron una altura de >5 y <10 m, 
seguido por el 26,17% de arboles que alcanzaron una de >10 y <15m. 23 arboles superaron los 40 m, 
que represento el 4,33% y el 30,69% de los arboles presento una altura de >15 y <40 m (Figura 2). 
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Figura 2 Distribucion de clases de altura de arboles presentes en parcelas (n = 531 individuos), Santo Domingo 

de los Colorados, Ecuador, 2007. 

De los 531 individuos inventariados las especies mas comunes fueron “Pechuga de gallina”, 
“Chfsparo” y “Pachaco”, y ninguna de las tres presento individuos en todas las clases diametricas. La 
“Pechuga de gallina”no presento individuos con dap > 50 cm, en cambio, “Chfsparo presento individuos 
con dap > 100 cm. “Pachacho” presento el 97,95% de sus individuos con dap >10 y < 40 cm. La especie 
“Chfsparo” no presento altura > 35 m, en cambio, “Pechuga de gallina” y “Pachaco” presentaron 
individuos con altura > 40 m (Figura 3). 





Figura 3 Distribucion de clases diametricas y clases de altura en las 3 especies mas abundantes en parcelas de 

bosque secundario y ripario (n = 531 individuos), 

Santo Domingo de los Colorados, Ecuador, 2007. 

La familia STERCULIACEAE presento el mayor promedio de altura con 25,96±10,32 m. La 
mayor area basal y volumen presento la familia SAPOTACEAE con 0,34±0,20 m 2 y 2,1 9±1 ,1 0 m 3 
respectivamente. La familia LECYTHIDACEAE obtuvo la menor altura, area basal y volumen con 
8,27±1 ,13 m, 0,02+0,00072 m 2 y 0,07+0,01 m 3 (Cuadro 5). 

3.4 Diferencias entre habitats 

En el bosque secundario se encontraron 282 arboles que corresponden a 43 especies distintas y 
en el bosque ripario se inventariaron 249 arboles que corresponden a 40 especies, agrupados en 23 
familias para los dos casos (Cuadro 5). 
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Cuadro 5 Habitat, numero de parcelas, arboles, especies y familias. 



Habitat 


N 9 de 

parcelas 


N 9 total de 

arboles 


N 9 total de 

especies 


N 9 total de familias 


Bosque 

secundario 


8 


282 


43 


23 


Bosque ripario 


8 


249 


40 


23 



No existieron diferencias estadfsticas en la riqueza de especies de las arboles presentes en los 
habitats, a pesar que la curva de acumulacion de especies (Figura 4) sugiere que los bosques 
secundarios son los mas diversos. En cada habitat existio un aumento progresivo de las especies 
conforme aumentaban los individuos inventariados. 




Figura 4 Curva de acumulacion de especies arboreas en habitats en base al numero total de arboles inventariados 
por cada habitat, Santo Domingo de los Colorados, Ecuador, 2007. 

La altura promedio, area basal y volumen fue mayor en el bosque secundario que en el bosque 
ripario pero no presentaron diferencias significativas (Cuadro 6). 

Cuadro 6 Promedio (±Se) de altura, area basal y volumen de arboles de acuerdo al habitat. 

Santo Domingo de los Colorados, Ecuador, 2007 
(L etras diferentes en la misma fila indican diferencias significativas, p < 0,05 ). 



Variable 


Bosque ripario 


Bosque secundario 


Altura total (m) 


14, 68+0, 57a 


1 5+0,69® 


Area basal (m 2 ) 


0,07±0,01a 


0,09+0,02® 


Volumen (m J ) 


0,62+0, 08a 


0,78+0,14® 



En bosque secundario 134 individuos presentaron un dap >10 y < 20 cm, categorfa en la que 
agrupo el mayor numero de arboles. 2 arboles alcanzaron un dap >90 cm, en bosque secundario y 146 
arboles presentaron un dap >20 y <90 cm. En bosque ripario 5 arboles presentaron un dap >70 cm, y 
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115 arboles alcanzaron un dap >0 y < 20 cm. Con un dap >20 y <70 cm, en bosque ripario se 
inventariaron 129 arboles (Figura 6). 
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Figura 6 Distribucion de clases diametricas de arboles presentes en parcelas de bosque secundario y ripario (n = 
531 individuos), Santo Domingo de los Colorados, Ecuador, 2007. 

En bosque ripario y secundario el 29% de los arboles inventariados tuvieron una altura > 5 m y < 
10 m. En bosque secundario y ripario el 52,83% y 61,45% de los arboles presentaron una altura > 10 m 
y < 40 m respectivamente. En bosque ripario el 2% de los arboles inventariados tuvieron una altura > 40 
m; en cambio el 6,38% de los arboles inventariados en bosque secundario tuvieron una altura > 40 m 
(Figura 7). 





□ Ripario 
■ Secundario 



CLASES DEALTURA 



Figura 7 Distribucion de clases por altura de arboles presentes en parcelas de bosque secundario y ripario (n = 531 
individuos), Santo Domingo de los Colorados, Ecuador, 2007. 

No se observaron diferencias significativas entre los habitats para los indices de Shannon (p = 
0,5825), Berger - Parker (p = 0.8979), Simpson (p = 0.8757), Margaleff (p = 0.6703) y Alpha (p = 
0.2542). Para Mackintosh (p = 0.0453) si se observo diferencia significativa entre los habitats (Cuadro 7). 

Estos resultados son similares a los encontrados en la zona de Rio Frio en donde no se 
reportaron diferencias significativas para los indices de Shannon y Simpson entre bosques densos 
intervenidos y bosques riparios (Chacon 2003). Sin embargo en otro estudio realizado en guanacaste, se 
determino diferencias estadisticas para los indices de Shannon, Simpson y Margaleff, entre los 
diferentes habitats, debido a la gran diferencia fisica existente entre habitats como cercas vivas y 
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bosques (Cardenas 2002) y en otro estudio desarrollado por Correa (2000), en Nicaragua se determine 
diferencias estadfsticas para el fndice de Alpha debido a la diferencia en la intervencion del hombre en 
los habitats en estudio como en bosque y cafe con bosque por lo que se pudo encontrar gran diversidad 
de especies en el bosque, pero en el habitat cafe con bosque pocas especies arboreas se encontraban 
en combinacion con el cultivo de cafe. 



Cuadro 7 Promedio indices (±Se) de diversidad y abundancia de especies de acuerdo al habitat. 
Santo Domingo de los Colorados, Ecuador, 2007 
Letras diferentes en la misma fila indican diferencias significativas, p < 0,05). 



indice 


Bosque ripario 


Bosque secundario 


CV% 


Shannon 


0,92±0,1 1 a 


1,00+0,07 a 


30 


Berger - Parker 


0,38+0,07 a 


0,36+0,07 a 


53,59 


Simpson 


0,23±0,06 a 


0,21+0,07 a 


82,22 


Margaleff M 


40,77±1 ,63 a 


39,64+2,02 a 


12,91 


Alpha 


5,63+1 ,36 a 


8,98+2,47 a 


77,09 


Mackintosh 


1 ,1 9+0,01 a 


1,15+0,01 b 


2,79 



4 CONCLUSIONES 

La diversidad de arboles de los bosques secundarios y bosques riparios de la zona de Santo 
Domingo de los Tzachilas es bastante alta, pero se puede estar perdiendo rapidamente por la extraccion 
de madera. 

Existe potencial para manejar los arboles existentes en los habitats, pero la regeneracion de mas 
arboles es muy baja. Por lo que es importante incentivas la proteccion de arboles en peligro de extincion 
y tratar de manejar los arboles en las tineas para tener beneficios para la conservacion y la produccion. 

5 IMPLICACIONES 

Se debe reconocer la importancia de los arboles en los habitats de las tineas para integrar estos 
elementos en los planes de conservacion y produccion. Se podrfa aumentar su valor para la 
conservacion: diversificando las especies presentes; aumentando su diversidad estructural y ubicandolas 
estrategicamente en el paisaje. Los arboles en las tineas deben y pueden ser utilizados tanto para fines 
productivos como de conservacion. 
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